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ABSTRACT 


The effects of radiation in GaAs solar cells has been extensively researched and 
the results of numerous investigation have yielded a considerable amount of information 
about the degradation in irradiated solar cells. This thesis establishes a novel method in 
which to use Silvaco's physically-based device simulator, ATLAS, to model the effects of 
radiation on solar cell output characteristics. A virtual model representing a single junc- 
tion GaAs solar cell was created in ATLAS. The effects of radiation were modeled using 
carrier trapping statements representing the defects associated with various fluence levels 
of 1 MeV electron radiation which were characterized with Deep Level Transient Spec- 
troscopy techniques. The resulting output characteristics of the virtual solar cell, ilumt- 
nated with a simulated AMO spectrum, were compared to published experimental meas- 
urements for GaAs solar cells of the same dimensions. The virtual solar cell demo n- 
strated a good correlation between the measured and virtual solar cell output characteris- 
tics and accurate representation of the spectral response. Complete ATLAS and 


MATLAB codes are included in appendices. 
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EXECUTIVE SUMMARY 


The military relies on satellite technology to fulfill communications and in infor- 
mational requirements. In current satellite configurations, the use of solar cells is vital to 
power production. The exposure of these solar cells to the degrading effects of radiation 
is of great concern in ensuring power requirements are fulfilled throughout the expected 
life of a satellite. There are a variety of methods currently utilized to predict the effects 
of radiation on solar cell output characteristics. Current methods of prediction concen- 
trate on determining a relations hip between the output characteristics observed over vari- 
ous radiation fluence levels for a given type of radiation. These methods require exhaus- 
tive testing over a variety of material parameters to determine the constants within the 
mathematical estimations for output characteristics. Thorough testing to establish mathe- 
matical relationships for materials such as Silicon and Gallium Arsenide has been ac- 
complished; however, in the evolution of solar cell design, new materials and solar cell 
configurations are continuing to be developed. The endeavor to mathematically model 
the output characteristics as a function of radiation type and fluence will continue to be 


an expensive effort due to the extensive fabrication and testing required. 


Previous theses by Andrew Bates and Max Green have established the ability of 
Silvaco's device simulator, ATLAS, to model solar cell configurations including mult+ 
junctions solar cells. This thesis is an extension of these works, as it determines a novel 
method to accurately predict the effects of radiation utilizing this physically-based device 
simulator. This method differs from the current methods of prediction due to the meth- 
ods of calculation employed by ATLAS. ATLAS uses a set of differential equations 
based on Maxwell's Laws to simulate the transport of carriers through a specified struc- 
ture. In order to develop a model to predict the effects of radiation, it was necessary to 
concentrate on the relationship between the ATLAS modeling parameters and the physi- 
cal results of radiation’s interaction with semiconductor materials. The permanent degra- 
dation of solar cells which have been exposed to radiation is the result of damage to the 
crystal structure of the semiconductor material. This damage can be measured in terms 


of traps, which are energy levels created within the bandgap of the semiconductor mate- 


XV 


rial. Through the use of Deep Level Transient Spectroscopy techniques, the characteriza- 
tion of traps induced by 1 MeV electron radiation in GaAs solar cell was investigated and 
published by Weng-Lyang Wang as part of his PhD dissertation at the University of Flor- 
ida in 1984. This trap information, along with additional trap parameters, published in 
1991 by Sheng S. Li and Robert Y. Loo in an article for Solar Cells, was used to model 
the effects of 1 MeV electron radiation in GaAs solar cells over various fluence le vels. 
The virtual solar cell parameters were obtained from an article published in JEEE Trans- 
actions on Electronic Devices, written by Robert Y. Loo, Sanjiv Kamath, and Sheng S. Li 


in 1990. 


ATLAS solar cell models were produced for four fluence levels of 1 MeV elec- 
tron radiation. The result of simulation demonstrated an accurate representation of the 
solar cell measurements published in 1990. There was, on average, less than 3% differ- 
ence between simulated and measured characteristics. The modeling of the virtual solar 


cell produced an accurate representation of the measured solar cell’s spectral response. 


XVI 


I. INTRODUCTION 


The military's use of satellites greatly enhances it's ability to communicate and 
gather information relative to the war fighter. Solar cells are today's primary source of 
power for satellite and space-based systems. These space-based systems rely on solar 
cells to meet overall system power requirements throughout their mission lives. Once in 
orbit, it is generally seen as prohibitively expensive, if not impossible, to repair or replace 
systems aboard spacecraft. Therefore, a considerable amount of research and planning 
are accomplished prior to launching systems bound for space. These systems are ex- 
pected to endure harsh launch and operating environments and ultimately accomplish 


their mission requirements [Ref. 1]. 


One of the primary environmental concerns affecting solar cells in space is expo- 
sure to radiation. The space radiation environment consists of cosmic rays, electroma g- 
netic radiation, Van Allen belt radiation, auroral particles, and solar flare particles [Ref. 
2]. Exposure to radiation over extended periods of time has been shown to degrade the 
power generation capability of solar cells [Ref. 2]. This reduction in supplied electronic 
power can adversely affect satellite missions, especially if the end-of-life power available 
does not meet mission power requirements, resulting in partial or complete mission fail 
ures. As the costs of space missions are very high, considerable effort has been and con- 


tinues to be expended to minimize the effects of radiation. 


Current methods used to predict the effects of radiation in solar cells revolve 
around the development of mathematical expressions that reflect the measured outputs in 
reference to radiation fluence levels [Ref. 3]. This process requires extensive testing to 
determine the relationship of the output to the fluence levels for a single energy level of 
radiation. This process is further complicated by the multitude of parameters that can be 


varied within a solar cell design. 


Green [Ref. 4] and Bates [Ref. 5] have experimented with modeling and optimiz- 
ing multi-junction solar cell designs using the ATLAS device simulation software by Sit 


vaco International [Ref. 6]. Their work demonstrated the usefulness of ATLAS as 


a design tool for solar cells. As an extension of these works, this thesis incorporates a 
novel approach to modeling and predicting the effects of radiation in solar cells designed 


in ATLAS. 


Chapter II provides a brief background discussion in semiconductor physics, solar 
cell operations and the nature of defects affecting solar cells. Chapter III describes the 
experimental measurements and values obtained from reference documentation. Chapter 
IV describes the software and modeling process utilized in simulating the referenced so- 
lar cells. Chapter V provides a comparison and discussion of the simulated and measured 
solar cell characteristics. Chapter VI notes the conclusions of this thesis and recommen- 
dations for future work in this area. Appendix A and B provide the ATLAS and 


MATLAB codes used to simulate the solar cells and radiation effects. 


Il. BACKGROUND 


This thesis presents a novel method to model the effects of radiation on solar 
cells. This chapter provides a brief discussion in basic semiconductor physics, solar cell 
operation and manufacturing processes. In addition, it also provides a brief description of 


radiation sources in space, their effect on solar cells, and the characterization of defects in 


semiconductors. 
A. SOLAR CELL BASICS 
if Semiconductor Physics 


There are several factors that govern the properties of semiconductors used in so- 
lar cells. Many of these properties are inherent to the atomic structure of the semiconduc- 
tor material. Silicon (Si) is one of the predominant semiconductors currently in use and 
will be used to elaborate the basic concepts of semiconductor physics. Si is a Group IV 
element which alludes to the four valence electrons of an isolated Si atom has in its outer 
shell [Ref. 7]. When in the vicinity of other Si atoms these valence electrons are shared 
to form covalent bonds between atoms in order to fill the atoms outer shells as shown in 


Fig. la. These bonds allow for the formation of orderly crystal structures. 
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Figure 1. Covalent bonds produced with (a) Group IV atoms, (b) Group III atoms, 
and Group V atoms [After Ref. 5.]. 


At a temperature of zero Kelvin, all valence electrons are contained in the valence 
band [Ref. 7]. When energy is imparted to these electrons, such as an increase in tem- 
perature, some of the electrons will break free from their parent atom and are free to 
move within the material. Free electrons will occupy the conduction band and the ab- 
sence of these electrons results in the formation of a hole in the valence band, represented 


by an empty circle in Fig. 1. In a pure, flawless material, there are no stable energy states 


3 


between the conduction and valence bands, this forbidden region is called the band gap. 
The energy required for an electron to transition from the valence band to the conduction 


band is referred to as the bandgap energy, E, [Ref. 7]. 


Conduction Band @ @ 





Valence Band @O @@ O 


Figure 2. Band gap diagram of a semiconductor [After Ref. 7.]. 

Controlling the number of electrons available for conduction is achieved through 
the addition of foreign atoms referred to as dopants. The selection of dopant atoms will 
determine if electrons are added or removed from the pool of electrons available for con- 
duction. Group III and Group V elements can be used to dope Group IV materials [Ref. 
5]. Atoms from Group III materials have three electrons in their outer shell and will form 
bonds with adjacent Si atoms. This will result in the formation of a hole as only seven of 
the eight available states in Si outer shell will be occupied as illustrated in Fig. 1b. At- 
oms from Group V materials have an extra electron which will be readily available for 
conduction as depicted in Fig. lc. Negatively doped Si, mtype, is the result of doping 
with Group V atoms while positively doped Si, p-type, is the result of doping with Group 
III atoms. It is also possible to produce materials with semiconductor properties by com- 
bining Group III and Group V elements such as InP or GaAs. Addtionally, alloys of dif- 
ferent mole fractions of Group III or Group V elements can be used to create ternary 


(three element) and quaternary (four element) semiconductors [Ref. 5]. 


A p-n junction is established when a p-type semiconductor is adjacent to an nh 
type semiconductor. The electrons in the mtype material will cross the junction to oc- 
cupy available positions in the p-type material and fill the atom’s outer shells as illus- 
trated in Fig. 3a. This transition of electrons across the junction will result in a positive 


region within the mtype and a negative region within the p-type material known as the 
4 


space charge region [Ref. 7]. The positive and negative nature of the space charge re- 
gion will create an electric field across the p-n junction [Ref. 7]. As this electric field in- 
creases, it will oppose the diffusion of electrons across the junction, stabilizing to the 


built-in potential, Eo, shown in Fig. 3b. 


Space Charge Region 








Neutral p-region Neutral n-region 


Electric Field ( E,) 


Group™ e- Group IV 
dopant dopant 


(a) (b) 


Figure 3. p-n junction depicting (a) the diffusion of electrons, and (b) the space 
charge region and built-in electric field [After Ref. 7.]. 

2. Solar Cell Operation 

The formation of the space charge region in a p-n junction and the associated 
electrical field are vital to the operation of solar cells. A solar cell in its simplest form is 
a p-n junction with contacts on both ends to extract current. As light enters the solar cell, 
photons will impart energy to the electrons within the space charge region [Ref. 9]. 
Those photons with energies greater than Fy will free an electron and create a hole, 
known as an electror-hole pair (EHP). EHP's created in the space charge region will be 
immediately separated by the electric field within the region. The hole will drift towards 
the negative end of the electric field while the electron will drift towards the positive end 
as shown in Fig. 4 (1). When the electron reaches the neutral n-type region it will make 
this region negative. Likewise, the hole will make the neutral p-type region positive 
[Ref. 7]. An open circuit voltage is then created between the positive p-side and the 
negative nside. EHP's created by photons outside the space charge region of the solar 
cell will diffuse based on carrier concentrations as there is no electric field with which to 
interact. The minority carrier of an EHP created within close proximity to the space 
charge region can travel the distance to the space charge region and be swept by the elec- 


tric field adding to the potential charge in the solar cell as shown in Fig. 4 (2 and 2'). 
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This distance, called the minority carrier diffusion length, is the mean distance an elec- 
tron, L., or hole, L;, can travel before recombining with an opposite charge [Ref. 7]. If an 
EHP is created outside the minority carrier diffusion length, the minority carrier will re- 
combine with an opposite charge and will not add to the accumulation of potential charge 
in the cell, Fig. 4 (3 and 3’). 
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Figure 4. Electron hole pair generation within the space charge region (1) ,minority 
carrier diffusion length (2 and 2’), and neutral regions of a solar cell(3 and 3’). 


Solar cell output characteristics can vary depending on cell configurations, types 
of material used, doping levels, and region thicknesses. In order to compare various 
cells, several standard parameters have been established. Many of these parameters can 
be expressed graphically on a current-voltage characteristics curve or /-V curve (see Fig. 
5). The limits of solar cell operations are set by the open circuit voltage, Voc, and the 
short circuit current, Jsc [Ref. 9]. The Voc is potential created by the accumulation of 
charge swept to either end of the space charge region while the solar cell is illuminated in 
an open circuit condition. The short circuit current is the measure of current expected in 
shorted condition. On the /-V curve, Voc is voltage corresponding to zero current and Isc 
is the current corresponding to zero voltage [Ref. 9]. The product of current and voltage 
pairs along the /-V curve correspond to an output power and the maximum output power, 


Pinax, 18 generally located near the knee in the /-V curve [Ref. 9]. 
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Figure 5. Typical J-V curve of a solar cell. 


Two other measurements by which solar cells are compared are fill factor, FF, 
and efficiency, ?. The FF is the ratio of P»,q, to the product of Isc and Voc [Ref. 9]. 
Graphically, the FF is a measure of the sharpness of the knee in the solar cells /-V curve. 
The efficiency of a solar cell is the ratio of P,,a, to the incident light power on the solar 


cell, Payait [Ref. 9]. The equations are [Ref. 9] 





FF = Fras (0.1) 
I Voc 
n=. (0.2) 


avail 


Another quantity that can be used to compare solar cells is the spectral response. 
The spectral response is the number of electrons generated at a given wavelength relative 
to the number of photons incident upon the solar cell at that wavelength, also known as 
quantum efficiency (QE) [Ref. 9]. The internal spectral response is defined as the num- 
ber of EHP's collected under short circuit conditions, /sc, compared to the number of pho- 
tons entering the cell [Ref. 9]. The external spectral response is the internal spectral re- 
sponse taking into account the reflection of light at the surface [Ref. 9]. Spectral re- 
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sponse can be regarded as a measure of how efficiently a solar cell converts a given 


wavelength of light to usable current. 
3. Growth Methods and Grown-in Defects 


In a search to produce higher efficiency solar cells, various methods of material 
growth have been developed. Each successive process has increased the purity and qual 
ity of the crystals produced. There are several methods for growing II-V compounds 
such as Gallium- Arsenide, GaAs. These can be separated into two basic categories: bulk 
growth and epitaxial growth. The two bulk growth methods include Liquid Encapsulated 
Czochralski (LEC) and the Bridgman method [Ref. 10]. The LEC grows a single crystal 
of GaAs from a molten liquid. The basic technique is modified due to Arsenic's volatility 
at Gallium's melting temperature. The liquid is covered with a layer of boron trioxide 
and a high pressure inert gas to prevent the evaporation of As. The Bridgman method 
uses a two-zone technique where a crucible is pulled from one side of a tube where the 
temperature is above the melting point of Ga to the other side where the temperature is 


reduced to accommodate As [Ref. 10]. 


Epitaxial growth methods are characterized by the growth of one crystal upon an- 
other [Ref. 10]. This growth can be of the same compound, known as homoepitaxy, or of 
a different compound from the substrate, heteroepitaxy. Various methods used in ept- 
taxial growth include Liquid Phase Epitaxy (LPE), Chemical Vapor Decomposition 
(CVD), Molecular Organic Chemical Vapor Deposition (MOCVD), and Molecular Beam 
Epitaxy (MBE) [Ref. 10]. 


Liquid Phase Epitaxy (LPE) deposits layers of material while they are in the liq- 
uid phase. When impurities are added to the liquid material, the melting point is reduced. 
As the resulting solution is cooled, new layers form on the seed crystal. As the crystal 
grows, the concentration of impurities in the liquid begins to increase, allowing a further 
drop in temperature. The rate of crystal growth is controlled by the rate of temperature 


reduction [Ref. 10]. 


Chemical Vapor Deposition (CVD), also known as Vapor Phase Epitaxy (VPE), 
grows crystals from gaseous compounds containing Ga and As. These compounds react 


on the substrate to form GaAs and byproducts of the remaining elements [Ref. 10]. 
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Metal Organic Chemical Vapor Deposition (MOCVD) is a specific type of CVD 
growth process. Metalorganic compounds are used to produce both the main elements 
and doping elements. This method is known to produce high quality material with large 


mobilities and minority carrier lifetimes [Ref. 10]. 


Molecular Beam Epitaxy (MBE) is a technique where gaseous beams of mole- 
cules interact with a crystal surface under ultra-high-pressure to form an epitaxial layer. 
The high standard of cleanliness and vacuum conditions, permit a very high level of pu- 


rity in the materials grown [Ref. 10]. 


The quest for high quality material is an attempt to reduce the number of defects 
in crystal structures and their effects on carrier lifetimes. Defects can be categorized as 
point, linear, or planar [Ref. 7]. Point defects may consist of interstitial atoms, vacancies 
and substitutional impurities in the crystal, Fig. 6a. Linear defects include edge and 
screw dislocations, Fig. 6b, while planar defects occur at grain boundaries, Fig. 6c. These 
defects are disruptions in the crystal pattern that permit energy states within the bandgap. 
These states are occupied by carriers during operation and affect carrier lifetimes, diffu- 
sion lengths, and increase the probability of recombination. Another source of defects 
within the solar cell crystal structure is the result of events external to the solar cell such 


as radiation. 





(a) (b) (c) 


Figure 6. Types of defects: (a) point, substitutional defect, (b) linear, dislocation de- 
fect, (c) planar, grain boundary defect [After Ref. 7.]. 


B. RADIATION 

iL Sources of Radiation 

Radiation of concern to solar cells in space can come from a variety of sources. 
These may include electrons and protons trapped in the geomagnetic field, radiation as- 
sociated with solar flares, and galactic cosmic ray radiation to a small extent [Ref. 3]. 
This radiation can be classified as one of two types: electromagnetic or particulate radia- 
tion. Electromagnetic radiation consists of light, X-rays, and gamma rays; particulate 
radiation consists of electrons, protons, neutrons, alpha particles, and others with a finite 
rest mass [Ref. 3]. Particle radiation can be characterized by the energy level of the par- 
ticle and fluence of the radiation encountered. The energy level is defined by units of 
electron volts (eV) which is equivalent to 1.6x10°'° joules. Fluence is defined by the total 
particles per area in a given time period with units of particles per cm’. 

2. Effects of Radiation on Solar Cells 

Radiation can interact with solar cells in a variety of ways but the primary damage 
mechanisms of interest in solar cell degradation are ionization and displacement damage 
[Ref. 3]. Ionization occurs when bound electrons are removed from their parent atoms. 
Radiation can affect solar cell outputs through several ionization related effects. These 
include the reduction of coverglass transmittance, darkening of transparent polymers, and 
altering of polymers mechanical properties [Ref. 3]. Ionizing radiation can also produce 
EHP in the same manner as light does but the energy level required to create an EHP is 
significantly greater than associated with light [Ref. 3]. The other mechanism, displace- 
ment damage, occurs when radiation dissipates energy to the solar cell, resulting in the 
displacement of an atom from its position within the crystal structure. When an atom is 
displaced, a vacancy is created, along with an interstitial atom, also known as a Poole- 
Frenkel defect, Fig. 7 [Ref. 11]. Higher energy radiation can have the effect of displacing 
multiple atoms in an elastic collision cascade [Ref. 11]. Radiation induced defects result 
in energy states within the band gap similar to those created by grown-in defects. These 
energy states result in trapping and recombination centers throughout the solar cell mate- 


rial [Ref. 11]. 
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Figure 7. Poole-Frenkel defect. 


3. Recombination and Capture of Electrons and Holes 

Figure 8 shows a series of electron transitions that will be used to illustrate re- 
combination and carrier capture. Fig. 8a represents the transitions associated with the 
absorption of light in the fundamental band and localized traps [Ref. 11]. Transition 1 
represents the generation of an electron hole pair as an electron is excited directly from 
the valence band. Transition 2 represents an electron excited from a bound state resulting 
in a free electron and a bound hole. Transition 3 is an electron excited from the valence 
band to a trap resulting in a free hole [Ref. 11]. Electrons and holes produced as a result 
of photo- ionization will exist until captured by a trap. Traps can exist as either trapping 
centers or recombination centers [Ref. 11]. Trapping centers exist at energy levels closer 
to the conduction or valance band. Carriers captured by trapping centers have a higher 
probability of being thermally excited back to a free state [Ref. 11]. Recombination cen- 
ters exist at energy levels closer to the center of the band gap and carriers captured have a 


higher probability of recombination with a carrier of the opposite charge [Ref. 11]. 
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Figure 8. Electron transitions: (a) photoionization, (b) carrier capture, and (c) re- 
combination [From Ref. 11.]. 


In Fig. 8b, the processes 4, 4', 5, and 5' represent the capture and thermal excita- 
tion of carriers in trapping centers [Ref. 11]. The transition 6 in Fig. 8b represents the 
capture of an electron by a recombination center and transition 7 is the capture of a hole. 
Temperature plays a significant role in whether a trap will act as a trapping or recombina- 
tion center. Figure 8c represents transitions associated with recombination. The transi 
tion 8 is the direct recombination of a free electron and a free hole. The release of a pho- 
ton is associated with this type of recombination but the probability of this event is low 
[Ref. 11]. The transitions 9 and 10 represent the more likely cases of capturing an elec- 
tron and hole, respectively [Ref. 11]. Excess energy from this transition can be dissipated 
through photons, phonons, or the Auger effect. The Auger effect occurs when the energy 
of two recombining carriers is transferred to a third carrier resulting in impact recombina- 


tion [Ref. 11]. 


C. CHARACTERIZATION OF TRAPS 

Traps can be characterized by various parameters to include the type of carrier 
they can be expected to capture, energy level, capture cross section, and the density of a 
specific trap encountered. 

1. Trap Type and Energy Level 

Traps are defined as electron or hole traps depending on which carrier they can 
expect to capture. Electron traps are expected to trap electrons from the conduction band 


and the trap energy level, E;, is measured from the conduction band to the local trap en- 
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ergy level. Hole traps are expected to capture holes from the valence band and their en- 


ergy level is measure from the valence band to the local trap level, Fig. 9 [Ref. 11]. 
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Figure 9. Carrier traps and energy levels. 


2 Capture Cross Section 

Carrier lifetime is not generally governed by direct recombination, but by the 
presence of defects [Ref. 11]. Recombination at traps is influenced by the occupation of 
these traps by carriers. The extent of occupation is related to the carrier density through 
the probability of capture by a vacant trap known as the capture cross section, s [Ref. 11]. 
The value of the capture cross section is determined by the potential near the trap. Traps 
whose electric charge attracts a carrier of opposite charge have large capture cross sec- 
tions, Fig. 10a. Neutral traps for a single displaced atom would have a capture cross sec- 
tion on the order of the atomic cross section of the displaced atom, Fig. 10b [Ref. 11]. 
Centers whose charge is the same as a carrier will be surrounded by a potential barrier 
causing electrostatic repulsion, resulting in a very low capture cross section as seen in 
Fig. 10c. A single displaced Arsenic atom within the crystal structure will be used to il- 
lustrate the nature of the capture cross section. Assuming the displaced atom creates a 
negative trap when displaced, the displaced atom will have a large potential for capturing 
holes. Once the trap has been occupied by a hole through the emission of an electron, it 
becomes neutral in charge and now has a smaller capture cross section for capturing elec- 


trons. 
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Figure 10. | Representation of capture cross sections; (a) opposite, (b) neutral, and (c) 
like charges with an attractive center [From Ref. 11.]. 


3. Trap Density 
The number of traps in a given volume is defined as the trap density, N,. There 


exists a linear relationship between trap densities and radiation fluence. As the fluence 
level increases, the density of traps also increases by a constant factor that is the product 
of atoms per unit volume, 7, capture cross section, s, and average displacements per 
primary displacement, ? [Ref. 3], given by 

N,=n,ov®. (0.3) 


As the number of particles irradiating the material increases, the quantity or den- 


sity of traps increases, but the type of damage created remains constant for a specific type 


of radiation at a given energy level. 


D. SUMMARY 
The background information introduced in this chapter is intended to provide the 


reader with the basic concepts and definitions to be used in the following chapters. The 
relationships and definitions provided are the basis for the development of a technique to 
predict the effects of radiation on a solar cell using ATLAS. This information also clari- 


fies the discussion of techniques and information annotated in the referenced expert 


ments. 
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Il. EXPERIMENTAL SOLAR CELL 


This chapter describes the physical structure and the measured the output charac- 
teristics of the referenced solar cells. It also presents the defects resulting from electron 
radiation and describes the experimental methods used to in characterizing the se defects. 
A. SOLAR CELL REFERENCES AND PARAMETERS 

The experimental data for this thesis comes from three published sources [Ref. 12, 
13, and 14]. The solar cells used for the experimental and structural data were described 
in Ref. 12, which is a study in radiation damage and annealing in GaAs solar cell's. Ref 
erence 12 provided information pertaining to the irradiation of the solar cells at room 
temperature with 1 MeV electrons at fluences ranging from 1x10'* to 1x10!° elec- 
trons/em’. The defect parameters were published in Ref. 13, which provided an overview 
of native and radiation induced defects in GaAs. Reference 13 is based in part on Ref. 
14, which is a study in grown-in and radiation induced defects and is the source of the 
experimental methods for defect characterization. 

1. Structural Parameters 

The solar cell from Ref. 12 was built as a p-on-n GaAs cell using the LPE growth 
method. The solar cell structure as shown in Fig. 11 consisted of three layers on an n'- 
GaAs substrate with a doping concentration of 5x10'’ cm*. The n-GaAs base was im- 
mediately above the substrate with a thickness of 10 um and had doping concentration of 
1x10'’ cm?. Above the base was a p-type GaAs emitter with a thickness of 0.5 um and a 
doping concentration of 1x10'* cm*. The window layer was made of AlGaAs and had a 
thickness of 0.5 um. The molar fraction of the AlkxGa;-xAs was listed as 0.95 and the 
doping concentration was 1x10'* cm®. The anode contact was made of AuZn with a 
thickness 300 nm and had a 4um layer of silver. The cathode contact was made of 
AuGeNi with silver overlay. The solar cell was produced with a Tantalum oxide, Ta2Os, 


antireflective coating but was not provided cover glass protection [Ref. 12]. 
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Figure 11. | Referenced GaAs solar cell structure. 


B. METHODS OF MEASUREMENT 

Many of the defects described in Ref. 13 were based on the work of Weng-Lyang 
Wang while at the University of Florida. Wang's PhD dissertation [Ref. 14] annotates 
several methods utilized to extract the damage parameters resulting from electron radia- 
tion. These included current-voltage characterization, thermally stimulated capacitance 
and deep level transient spectroscopy. 

1. Current- Voltage (I-V) Characterization 

Measuring the current-voltage characteristics of a p-n junction diode, in the for- 
ward bias region, can reveal information concerning the recombination process within the 
space charge region [Ref. 14]. The total current in a p-n junction diode is a combination 
of the diffusion current and the generation-recombination current in the space charge re- 
gion. When the generation recombination current dominates the current expression the 
effective carrier lifetimes can be derived from the ma gnitude of the generation- 
recombination current. The density of traps, N;, and capture cross section, sn, were calc u- 
lated using the carrier lifetimes defined with reference to the average thermal velocity, 


ih, aS 


Tt, =———_.. (3.1) 


2. Thermally Stimulated Capacitance (TSCAP) 


Another technique used to calculate trap densities is known as the thermally 
stimulated capacitance (TSCAP) measurement [Ref. 14]. TSCAP was accomplished by 
reverse biasing a p-n diode and then cooling the temperature to that of liquid nitrogen (77 
K). The diode was briefly zero biased and then reversed biased [Ref. 14]. This allowed 
the majority carrier traps to be filled. The temperature was then raised to around 400 K. 
A thermal scan of the capacitance, C, verses temperature, 7, was plotted. If carrier emis- 
sion were taking place a capacitance step was observed on the C vs. T plot. The ampli- 


tude of the capacitance step is proportional to the trap density. 
as Deep Level Transient Spectroscopy (DLTS) 


Deep Level Transient Spectroscopy (DLTS) was introduced by Lang in 1974 as a 
technique to study deep levels in semiconductors [Ref. 3]. Deep levels in the depktion 
region of a p-n junction can be filled with charge carriers by an electrical pulse [Ref. 14]. 
The capacitance transient is monitored as the temperature is increased at a constant rate 
and a DLTS spectrum is produced. This spectrum is used to extract trap parameters 
within the tested region. These characteristics include the trap energy, E;, capture cross 


section, $ y,», and trap density, N;. An example of the DLTS spectrum is shown in Fig 12. 
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Figure 12. | DLTS scans of electron traps for two LPE GaAs samples [From Ref. 14]. 
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Cc. EXPERIMENTAL MEASUREMENTS 

L Measured Trap Values 

A thorough description of electron and hole trap information for both native and 
radiation induced defects was collected by S. S. Li and R. Y. Loo [Ref. 13]. Trap data 
was available for various fluence levels and temperatures. The values listed in Table | 
are traps measured in specimens irradiated at 25°C across three fluence levels of electron 
radiation [Ref. 12]. The range of data at this temperature provided enough detail to accu- 
rately characterize the effects of electron radiation over a wide range of fluences. Table 1 
is partitioned into electron and hole traps with respect to electron fluence levels. Each 
trap is identified by energy level, density, and capture cross section for the specified car- 
rier. Electron traps were measured in the n-GaAs region while the hole traps were meas- 
ured in the p-GaAs regions of the mesa diodes. Data marked with an asterisk were not 
provided in the references and was subsequently estimated based on traps with similar 
energy levels or the linear relationships between fluence and density, as previously dis- 


cussed. 


Table 1. Experimentally measured electron and hole traps in n and p-GaAs [From 
Ref. 13]. 


Fluence | Electron es in n-GaAs Hole traps in p-GaAs 


F, E, E; Nt Sp 
= (eV) pe (em?) (eV) (cm*) (cm?) 


eel [82x07 [2eio—[ 
| Ec-0.71 [70x10 | 83x10" | | 
| Ec-0.90 | 88x10" | 30x10" | | 
arse are as er es a 
Ec-0.At 
eeEeEe 
a A PE PTC 
ee Sneren ere kc Sea 
7.2x10. 
Ee or eo 2.2x10 


* interpreted values based on similar trap energy levels and linear relationships 
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Figure 13 is provided to confirm the linear relationship of trap densities with flu- 
ence levels. The linear relationship holds for those traps with data available over three 
fluence levels. The exception to this would be the hole trap at an energy level of 0.71 eV 
above the conduction band at a fluence level of 10'° e/cm’. This trap density appears to 
be higher than would be expected in relation to the trap measurements at the two adjacent 


fluence levels. 
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Figure 13. | Experimentally measured trap densities in relation to fluence levels [After 
data from Ref. 13.]. 

Za Measured Output Chracteristics 

The data available for solar cell output characteristics was normalized to the pre- 
radiation output values. Figure 14 is a plot of measured Isc, Voc and the fill factor for 
fluence levels up to 1x10'° e/em’. The output characteristics are shown to degrade with 
an increase in radiation fluence. The Jsc displays a 40% reduction in current from its pre- 
radiation value while Voc displays approximately a 15% drop in voltage [Ref. 12]. The 
fill factor shows little degradation over the range of fluence values. Graphically, this 
would be seen as a decrease in the size of the /-V curve while the general shape of the 


curve remains constant. 
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Figure 14. Normalized short circuit current, open circuit voltage, and fill factor [After 
experimental data from Ref. 12.]. 


The maximum power and solar cell efficiency are shown in Fig. 15. The Pina, 
prior to irradiation was 21.25 mW and drops to 53% of it original value at a fluence level 
of 1x10'° e/em? [Ref. 12]. The solar cell efficiency has a pre- irradiation value of 15.7 % 


and deteriorates to 8.3 % efficiency at the highest level of fluence [Ref. 12]. 
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Figure 15. Experimentally measured P,,,, and efficiency in relation to radiation flu- 
ence levels [after data from Ref. 12]. 
D. SUMMARY 
The output characteristics and trap parameters provided by the referenced docu- 

ments provide a basis for comparison of ATLAS simulations in the following chapter. 
The bulk of the information gleaned from the references is in adequate detail to produce a 
comparable model in ALTAS. The experiments were specifically tailored to support the 


subjects of the referenced articles. A lack of detail is provided in the range of error for 
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measurements and number tested. A more thorough investigation into differences in 
simulated and measured values would be possible had the data been specifically tailored 
for the purpose of validating the ATLAS model. The following chapter provides a brief 
description of ATLAS and a detailed discussion of modeling an equivalent GaAs solar 


cell. 
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IV. SOLAR CELL SIMULATION 


This chapter presents the radiation modeling approach and discussion of parame- 
ter selection decisions. The measured solar cell parameters and output characteristics dis- 
cussed in the previous chapter were the benchmarks for the simulations. The measured 
trap details provided the foundation for modeling radiation damage in the virtual solar 
cells. 

A. THE SILVACO ATLAS DEVICE SIMULATOR 

ATLAS is a physically based two- and three-dimensional semiconductor device 
simulator [Ref. 6]. Physically based device simulators predict the electrical behavior as- 
sociated with the physical structure under specified biased conditions. ATLAS achieves 
this through partitioning the specified device into a two or three-dimensional mesh grid. 
ATLAS then applies a set of differential equations, based on Maxwell's Laws, to the 


mesh to simulate the transport of carriers through the structure [Ref. 6]. 


Interface with ATLAS is provided through Silvaco's Deckbuild run-time envi- 
ronment via an ASCII text file. Structure and command inputs to ATLAS are required to 
simulate a device [Ref. 6]. The structure data contains the dimensions, material charac- 
teristics, and mesh grid parameters of the device desired. The command file identifies the 
commands ATLAS will execute in order to provide the desired output. The order in 
which statements are specified to ATLAS is significant and must conform to the follow- 
ing precedence: structural specifications, material models specifications, numerical 
method selection, solution specification, and then results analysis [Ref. 6]. Additionally, 
the order within each of these specifications is also significant. For example, when iden- 
tifying the structure specifications, the required command order is mesh, region, elec- 
trodes, and then doping. Failure to adhere to the specified order will result in ignored 
statements, incorrect simulation results and/or program termination. The ATLAS code 
for simulating the various models used for this thesis is provided in Appendix A. 

B. THE VITUAL SOLAR CELL MODEL 

The experimental solar cell’s material and dimensions were modeled as specified 

in Chapter III with a few exceptions. Some parameters, not identified by the references, 


were estimated while other parameters could not successfully be modeled on the current 
23 


version of ATLAS. A separate solar cell model was created to represent each level of 
radiation. The damage as a result of radiation was represented through a series of trap 
statements within the material model specification portion of the code. 

1. ATLAS Calculations 

ATLAS uses a mathematical model that operates on any semiconductor device. 
This model consists of a set of fundamental equations which link together the electro- 
static potential and carrier densities [Ref. 6]. The primary equations derived from Max- 
well’s laws consist of Poisson’s Equation, the continuity equations and transport equa- 
tions. Poisson’s Equation relates variations in electrostatics potential to local charge den- 
sities 


divie Vy)=-p, (4.1) 


where ? is the electrostatic potential, e is the local permittivity, and ? is the local space 


charge density. The electric field is obtained from the gradient of the potential 
E=Vy, (4.2) 


where E is the electric field. The continuity and transport equations describe the way 
carrier densities are calculated as a result of transport, generation and recombination 


processes [Ref. 6]. The continuity equations for electrons and holes are 


m Hair eG: -R,, (4.3) 
eG 

op 1. = 

a = ai, +G,-R,, (4.4) 


where n and p are the electron and hole concentrations, J, and J, are the current densi- 
ties, G, and G, are the generation rates, K, and R, are the recombination rates, and q is 
the electron charge. 

2. Virtual Solar Cell Parameters 

The virtual solar cell shown in Fig. 16, was modeled, when possible, according to 
the dimension and material characteristics specified in Chapter III. The reference docu- 
mentation does not identify the width of the anode contact or the thickness of the cells 


substrate. When creating the virtual solar cell structure, the anode contact width was 
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modeled to be 8% of the cell's surface area and the substrate thickness was modeled as 
200 um. Additionally, the virtual cell was designed to produce an output for a 1 cn cell 


area, where the referenced cells had areas of 4 cm. 





Figure 16. — Structure of the virtual solar cell produce in Tony Plot. 


The referenced solar cells were illumination by an unspecified Air Mass Zero 

(AMO) spectrum simulator. The virtual solar cell was illuminated with the most current 
AMO spectrum file, Fig. 17, identified by the Renewable Resource Data Center, National 
Renewable Energy Laboratory (NREL) [Ref. 15]. The virtual solar cell contacts were 
modeled as pure gold for optical characteristics only; no resistance was associated with 
the modeled contacts. External lump series resistance simulations were unsuccessful on 
the current version of ATLAS. A five percent reflection of incident light upon the first 
layer was used to model the typical effects of the solar cell's Ta2Os antireflective coating 


[Ref. 5]. 
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Figure 17. Air Mass Zero spectrum file used in model illumination [From Ref. 15.]. 


Significant interaction occurs within the space charge region and the minority car- 
rier diffusion length. To ensure adequate calculations over these regions, it was neces- 
sary to build a solar cell structure that specifically identified the space charge region. The 
static space charge region was identified in the model based on the calculated depletion 


region widths from [Ref. 7] 


w= A (4.5) 


GN Np 


Ee | 
Dl- 
. 


where W is the depletion width, e is the local permittivity, N4 is the acceptor doping con- 
centration, Np is the donor doping concentration, EF, is the built in potential, and g is the 
charge of an electron [Ref. 7]. The complete depletion width is the total sum of the p- 


and mregion depletion region widths, shown by 
W=W,+W,, (4.6) 


where W,, is the width of depletion region in the nregion and W, is the width of the de- 


pletion region into the p-region. These widths were modeled to reflect the depletion 
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depths into the n and p-regions from the charge on one side of the depletion layer given 


by 
|Q|=N,W,A=4N.W,A, (4.7) 


where A is the area of the junction [Ref. 7]. Manipulation of Eqns 4.6 and 4.7 can be 
used to obtain the depletion widths corresponding to the p-GaAs emitter and n GaAs base 


by 


Ww 
W,=———_~, (4.8) 
N D 
1 + emitter. 
ND see 
Ww 
WwW =—~ (4.9) 
N 


where Np is the doping concentration of the specified region. The space charge regions 
identified by these calculations ensure an adequate number of mesh nodes are identified 
to provide the accurate calculation of carrier transport for this vital region. Figure 18 de- 
picts close-up views of the mesh within the p-n junction of the virtual solar cell for both 
identified, Fig. 18a, and unidentified, Fig. 18b, space charge regions. In Fig. 18a, region 
1 is the contact, regions 2 and 3 are vacuum, and region 4 is the AlGaAs window. The 
emitter is divided into two regions, region 5 is the neutral p-GaAs region of the emitter 
and region 6 is the space charge region of the emitter. Similarly, the base is divided into 
the base space charge region, region 7, and the neutral n-GaAs base region, region 8. The 
regions in Fig. 18b are similarly identified but without the emitter and base space charge 
regions. A comparison of the mesh layouts reveals a mesh with a significantly higher 
node density in the space charge region of the solar cell with an identified space charge 
region. In Fig. 18b, it can be seen that the area corresponding to the space charge region 
in the emitter, Fig. 18a (6), is completely omitted in the mesh layout and as a result, the 


vital region would be ignored during carrier generation and recombination calculations. 


Zs 





(a) (b) 


Figure 18. | Mesh layout in the vicinity of the p-n junction in the virtual GaAs solar 
cell model with (a) specified and (b) unspecified space charge regions. 

x Material Models 

ATLAS permits the user to identify a variety of models [Ref. 6] for calculating 
carrier mobility and recombination. The Concentration Dependent Low Field Mobility 
model (CONMOB) was used to model the doping-dependent low-field mobilities of elec- 
trons and holes in GaAs at 300K. Through CONMOB, ATLAS provides empirical mo- 
bility data relative to doping concentrations specifically tailored for Si and GaAs [Ref. 6]. 
The recombination models utilized were the optical (OPTR) and the Shockley-Read-Hall 
(SRH) recombination models [Ref. 6]. The OPTR model was designed to be used with 
direct band gap materials such as GaAs. SRH recombination model takes into account 
the electrons being emitted or captured by donor and acceptor- like traps. ATLAS also 
provides a SRH Concentration- Dependent Lifetime (CONSRH) model where carrier life- 
times become a function of impurity concentrations [Ref. 6]. The application and limita- 
tions of these SRH models are discussed in the Analysis of Results Chapter of this thesis. 

4. Trap Implementation 

When implementing traps and defects ATLAS utilizes Poisson’s equations in the 


form 


dive Vy) =q(n—p-N, + N,)+@Q,, (4.10) 
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where n and p are the carrier concentrations, Np and Ny are the donor and acceptor con- 
centrations, and Qr is total charge caused by the presence of traps. Q7, is added to the 


right side of Poisson’s equation [Ref. 6] with the value defined as 


QO, = q(p, —n,), (4.11) 


where n; and p; are the densities of the trapped charge. These trapped charge densities are 


defined by, 
n,=N,F (4.12) 

and 
D, =N,F (4.13) 


where WN, is the density of traps, and F’, and F;, are the probabilities of occupation by carri- 
ers. The probabilities of occupation are based on the capture cross section and emission 


rates of the individually defined traps. 


Based on Eqn. (4.10), ATLAS is capable of simulating a perfect crystal structure 
by setting Q7 equal to zero. However, this condition does not physically exist in large 
structures such as a solar cell. Impurities and defects associated with crystal growth are 
always present. To account for these defects, trap statements were implemented in the 
pre-irradiated solar cell model. Trap values prior to the application of radiation were not 
known for the referenced solar cells. Based on the LPE growth process, a trap energy 
level of 0.71 eV was selected for the pre-radiation model [Ref. 13]. The densities of 
these traps was interpreted based on the linear relationship between trap density and flu- 
ence level. The curve fit between data at fluence levels of 10'* and 10’° provides a linear 


relation ship in the form of 
N, =0.027®, +3x10"', (4.14) 


where WN; is the trap density and ? . is the electron fluence level. From Eqn. 4.14, the trap 
density for the pre-radiation model was determined to be of 3x10'! /em’. Traps in all 
other models were implemented in accordance with the referenced findings and previ- 


ously annotated interpolations. 
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Trap statements in ATLAS have a specific format and required information fields 


[Ref. 6]. A typical trap statement is expressed in the following form: 








trap acceptor e.level=0.71 density=le6 sign=4.0e13 sigp=4.0e15 region=5 


The trap statement is initiated by the term trap. Traps are identified as acceptor- or do- 
nor-type traps. ATLAS defines acceptor traps as "hole traps" measured from the conduc- 
tion band to the hole trap [Ref. 6]. Donor traps are defined as "electron traps" measured 
from the valence band to the electron trap, as shown in Fig. 19. This reveals that ATLAS 
utilizes holes in the conduction band and electrons in the valence band which is opposite 
to the definitions typ ically used by physicists [Ref. 11]. However, this program is meant 
to be integrated with complex electrical circuits and hole flow is in line with electrical 
engineering practices [Ref. 6]. As such, hole traps are implemented in the p-GaAs re- 
gions as acceptor traps while electron traps are implemented in the n-GaAs regions as 
donor traps. 

Conduction Band 


Eo 
Acceptor Trap 


Intrinsic Fermi Level 


E, 
Donor Trap cis 
E 


Vv Valence Band 


Figure 19. Representation of Silvaco trap definitions [From Ref. 7.]. 


Capture cross sections are labeled as sign or sigp referring to the charge dis- 
played by the trap. Prior to a trap capturing an electron or a hole, the capture cross sec- 
tion of a trap was implemented as that value measured experimentally for an oppositely- 
charged carrier [Ref. 6]. Once the trap has captured a carrier, it now becomes neutral and 
will have a smaller capture cross section [Ref. 6]. Measurements of neutral trap capture 
cross sections were not available in the referenced documentation. Each trap may be 
made up of a variety of situations. Some combinations may include a vacancy, an inter- 


stitial atom, multiple interstitials, multiple vacancies, or a combination of vacancies and 


30 


impurities [Ref. 14]. As a result, it is not known what exactly is causing a specific trap 
induced by radiation. It was speculated that the capture cross section of a neutral trap 
would be some constant ratio of the measured capture cross section value. This ratio was 
determined to be 1x10 through iterative simulations. To further elaborate, hole traps 
would display a greater capture cross section when negatively charged, represented in the 
value of the parameter sign, while having a smaller capture cross section when neutral 
represented in the parameter sigp. The proper implementation of traps was confusing to 
decipher with regards to the integrating of physics based data with an electrical engineer- 
ing-based simulation. To clarify the implementation of traps, an example using charac- 
teristics one would encounter reading physics documentation and implementation of this 


data in ATLAS code, is provided. 


A hole trap measured in the p-GaAs region with an energy level of E;, a measured 
hole capture cross section of sp, and a measured trap density of N; would be implemented 


in the following manner: 





trap acceptor e.level=F, density=N, sign=s, sigp=(s,) *le-2 region= 
Similarly, an electron trap measured in the n-GaAs region with an energy level of E;, a 
measured electron capture cross section of s ,, and a measured trap density of N; would be 
implemented in the following manner: 


trap donor .level=E, density=™, sign=(s,)*le-2 sigp=s, region= 





2: Results of Simulation 

The results of simulation are shown in Table 2 with reference to electron fluence 
levels, ? -. Short circuit currents and open circuit voltages were extracted from ATLAS 
output data. Maximum power, fill factor, and efficiency were calculated based on stan- 
dard definitions. The virtual solar cell modeled with ATLAS demonstrates a 15% drop in 
open circuit voltage and 35% drop in short circuit current at the highest fluence levels. 
The model shows a 50% loss in efficiency and maximum power with a small degradation 


in fill factor. 
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Table 2. Output characteristics of the virtual GaAs solar cell implemented in 
ATLAS. 


Pre-Rad ? .=10" ? .=10 ? .=10" 
e/em” e/em e/em” 


26.44 24.73 21,31 17.13 





The degradation in the /-V curve as a result of electron radiation is shown in Fig. 
20. The consistent shape of /-V curve was anticipated and is an indication of the fill fac- 
tor's consistency over changes in fluence levels. The maximum power points are anno- 
tated with diamonds on their respective curves. These results are consistent with typical 


experimental results available in literature [Ref. 16]. 


Prior to irradiation 





Current, mA 





0.0 0.2 0.4 0.6 0.8 1.0 1.2 
Voltage, V 


Figure 20. | Current- Voltage curve for the virtual GaAs solar cell as function of elec- 
tron fluence level. 


32 


The spectral response curves of the virtual solar cell are sho wn in Fig. 21. The 
pre-irradiation curve shows a peak value of 81% at a wavelength of 0.65 um. The dete- 


rioration of the spectral response is shown with reference to the associated fluence levels. 
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Figure 21. | Spectral response of the virtual solar cell with respect to radiation fluence 
levels. 


Figure 22 shows a series of the virtual solar cell model structures overlaid with 
the photogeneration rates corresponding to the annotated wavelength. The EHP genera- 
tions at shorter wavelengths occur uniformly across the cell. As the wavelength in- 
creases, the photogeneration increases deeper in to the solar cell. The shadow of the gold 
contact can be seen in the middle of the cell structure and it can be seen that gold at these 
small thicknesses permits the transmission of some light into the cell. The numbers listed 


in the are the Logio of the photogeneration rates for the given wavelength 


a2 





Figure 22. Illumination of the virtual solar cell model and the resulting photogenera- 
tion rates with respect to various wave lengths. 


The generation of photocurrent is only part of the output associated with an illu- 
minated virtual solar cell model. Figure 23 shows the recombination rates associated 
with annotated wavelengths. At shorter wavelengths the highest recombination occurs in 
the AlGaAs window layer. As the wavelength increases, the recombination rate increases 
deeper into the solar cell. One can also notice that the space charge region has a conss- 


tently low recombination rate due to the inherent electrical field in this region. 
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Figure 23. —_ [llumination of the virtual solar cell model and the resulting recombination 
rates with respect to various wavelengths. 


C. SUMMARY 

The parameters and techniques used to model the virtual solar cells in ATLAS 
were presented and justified in this chapter. The output characteristics of the virtual solar 
cells were shown to degrade with the implementation of traps associated with the three 
levels of electron fluence. The following chapter compares these results to the measured 
values of the referenced documentation and expounds upon possible causes of the more 


significant deviations in values. 


35 


THIS PAGE INTENTIONALLY LEFT BLANK 


36 


V. ANALYSIS OF RESULTS 


In this chapter the results of the previous chapter are compared to the measured 
results presented in Chapter III. Some of the differences in output characteristics are dis- 
cussed and potential causes are identified. 

A. DIRECT COMPARISON OF OUTPUT CHARACTERISTICS 

The simulation results of the virtual solar cell model are tabulated and compared 
to the experimental values in Table 3. There is, on average, less than 3% error between 
measured and calculated parameters. Simulated and measured values of Isc, Voc, and fill 
factor are compared in Fig. 24 while maximum power and efficiency are compared in 
Fig. 25. In both graphs, the simulated data is represented by discrete points, while the 
measured data is represented by lines derived from several points selected from the origi- 
nally published curves [Ref. 12]. The fill factor, as expected, is consistently high as a 
result of modeling with no resistance in the contacts. 


Table 3. Comparison of virtual model and experimentally measured output charac- 
teristics. 


Parameter 1 Mev Fe = 10'e/cem 


% [Ref. 12] % 
| = | 094 | 0.936 | 0.48 | 
| = | 0.96 | 0.973 | 130 | 


1 Mev Fe =10°e/cem 1 Mev Fe = 10'°e/em 
Normalized Values Normalized Values 


ai 
Meas. | Sim Diff, | Meas. Sim Diff, 
[Ref. 12] % [Ref. 12] % 
0.806 | 1.70 | 0.60 | 0.648 | 8.00 | 
; 1 





1 
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Normalized parameter 


Figure 24. 
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Figure 25. 
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Comparison of the virtual and experimentally measured solar cell's Pinax 
and efficiency with respect to radiation fluence levels. 


The results of the virtual solar cell for the fluence level of 10!° e/cm’ demonstrate 


the greatest disagreement in maximum power. Both /sc and Voc are 2% below the meas- 


ured values and the maximum power is 8% below the experimental measurement. This 


difference could be due to the inaccuracy of trap density information for this fluence 


level. As stated in Chapter II, the plot of trap densities showed a higher than expected 


value for the trap density at the hole trap with an energy of 0.71 eV. This energy level 


corresponds to a recombination center and would account for lower output characteristics 


at this fluence level. Reference sources provided no indication to the range of error in the 


measurements for trap densities. 
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At a fluence level of 10'° e/cm” there is an 8% difference in measured and simu- 
lated short circuit current. It is speculated that a programming limitation in the impk- 
mentation of recombination models may be the source of this difference. GaAs has dem- 
onstrated trap induced carrier lifetime variations that are dependent upon doping concen- 
trations [Ref. 17]. Shallow traps, near conduction and valence bands, demonstrate a re- 
duction in carrier lifetime as the ratio of the doping concentration to the trap density is 
decreased. Deep level traps, near the intrinsic Fermi level, have been shown to be inde- 
pendent of doping concentrations. ATLAS has the ability to account for a concentration- 
dependent recombination model with the CONSRH model statement. However, recom- 
bination models are identified over specific regions and cannot be applied to individual 
trap statements within a region. Simulations with the CONSRH model produced a short 
circuit current value that was 9% below the measured output value at a fluence level of 
10'° e/em’ while the SRH model was shown to produce a value 8% above the measured 
output value. A more accurate calculation may be obtainable if recombination models 
could be applied to specific trap statements. 

B. SPECTRAL RESPONSE 

Figure 26 is a graphical representation of the experimentally measured spectral 
response of an actual GaAs solar cell with reference to electron fluence for wavelengths 
of light ranging from 0.4 um to 0.9 um. These values are normalized with respect to the 
maximum pre-radiation spectral response value. Similarly, the virtual solar cell's spectral 
response for the same radiation fluence levels are shown in Fig. 27. Both experimental 
and simulated spectral response curves show peak values around a wavelength of 0.65 
um. The peak values on both plots are approximately 98, 88 and 72 percent of the origi- 
nal response for fluence levels of 10'*, 10’° and 10'°e/cm’, respectively. Significant de- 
terioration in spectral response can be seen at the higher wavelengths where both curves 
show similar reductions in the longer wavelengths on the curves. Simulations in ATLAS 
produce an accurate representation of measured spectral response with respect to radia- 


tion fluence levels. 
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Figure 26. | Measured spectral response as a function of fluence [From Ref. 12.]. 


100 








Spectral Response, % 
Nn 
So 








0.4 0.5 0.6 0.7 0.8 0.9 
Wavelength, Um 


Figure 27. | Simulated spectral response as a function of fluence. 
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C. SUMMARY 

This chapter presented and compared several experimentally measured output 
characteristics with corresponding virtual solar cell model characteristics. A method to 
accurately predict the effects of radiation, based on the trap levels detected in the material 
after electron irradiation, was demonstrated. The following chapter presents the conclu- 


sions and recommendations for future study based on this work. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


The radiation modeling methodology demonstrated in this thesis utilizes measur- 
able trap data to predict the degradation in output characteristics of solar cells. This dif- 
fers from current methods of modeling and prediction which rely on extensive testing and 
data accumulation to develop empirical mathematical relationships between radiation flu- 
ence levels and solar cell output characteristics. The demonstrated method may reduce 
the testing required to characterize radiation damage in solar cells by capitalizing on the 


linear relationship between radiation fluence levels and trap density. 


This thesis used published experimental information and measurements to virtu- 
ally model the dimensions, materials and the effects of radiation on a p-on-n GaAs solar 
cell. The experimental measurements for this thesis were obtained from three primary 
publications. The first [Ref. 12] source provided the measurements for the output charac- 
teristics and the degrading effects of electron radiation. The remaining sources [Ref. 13 
and 14] provided detailed trap characterization for grown-in and radiation induced de- 
fects in GaAs. Using this information a virtual solar cell model was created using Sit 
vaco's ATLAS a physically based device simulator. Virtual GaAs solar cell models were 
created based on the documented dimensions, material and doping concentrations. The 
published trap information was used to model the effects of radiation in the virtual solar 
cells. Simulations with the virtual soar cell model illuminated with AMO demonstrated 
an accurate representation of measured solar cell current, voltage, and spectral response 
values. The virtual solar cell demonstrated an accurate degradation in output characteris- 
tics when traps associated with corresponding fluence levels were implemented. The vir- 
tual solar cell model demonstrated accurate degradation over four fluence levels of elec- 


tron radiation. 


This method of modeling the effects of radiation can be used to predict the effects 
of radiation that produces uniform damage throughout the material [Ref. 19]. Electrons, 
neutrons, and high energy protons can be modeled using uniform trap implementation. 
Low energy proton radiation is a notable exception as it does not introduce damage unt 


formly within a semiconductor material [Ref. 12]. The penetration depth of low energy 
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protons is highly dependent upon energy levels resulting in a very complex trap density 


profile for which information was not immediately available. 


This methodology may also prove useful in modeling the effects of annealing in 
solar cells. Annealing occurs when energy is added back to the crystal either through 
temperature or current injection. This energy allows for some displaced atoms to return 
to their proper placement in the crystal and thereby reducing the density of traps. Similar 
to radiation fluence, there is a linear relationship between the reduction in trap density 
and the time a solar cell is held at constant temperature [Ref. 12]. This relationship may 
prove useful in predicting solar cell annealing behavior with minimal testing require- 


ments [Ref. 20]. 


This methodology shows the greatest potential with regards to multt junction so- 
lar cells. ATLAS may provide the ability to predict the effects of radiation over the en- 
tire structure of a multtjunction cell based on information obtained from individual di- 
odes of cells. The trap information required to model the effects of radiation can be ob- 
tained by irradiating separate p-n junctions produced of the same materials that make-up 
a multtjunction cell. The resulting trap information can be used to model the effects of 
radiation in the individual junctions and predict the degradation in the overall solar cell. 
Additionally, when used in conjunction with published optimization routines [Ref. 5], 
this methodology may be used to optimize mult junction solar cells for power based on 


the projected end-of-life radiation exposure. 


Physically based device simulation demonstrates promising potential to aid in the 
advancement of solar cell design [Ref. 20]. It can reduce the cost of design and produc- 
tion by minimizing testing requirements and optimizing structure design. This method of 
simulation may also provide a vehicle to exploit previously unused parameter relation- 
ship that may prove both accurate and insightful. Ultimately, the ability to reduce costs 
in design, testing and fabrication may aid in the proliferation of solar cells as a widely 
used renewable energy source. This work has been accepted for presentation at the 31st 


IEEE Photovoltaic Specialist Conference in Orlando, FL [Ref. 19 and 20]. 
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APPENDIX A. ATLAS INPUT DECKS 


This appendix provides the ATLAS input decks used for this thesis. Section 1 
provides the input deck for the pre- irradiation solar cell model. Section 2 provides the 
trap statements for the remaining fluence levels. 

1. PRE-IRRADIATION SOLAR CELL INPUT DECK 

The following input deck was used to simulate the pre- irradiation GaAs solar cell 
model. The code provides for the calculation and plotting of the current-voltage charac- 
teristics curve, the dark current-voltage characteristics curve and extraction of data to plot 
the spectral response curve. 


go atlas 

mesh width=200000 

#X—-Mesh 

x.mesh loc=-250 spac=62.5 
x.mesh loc=-20 spac=2.5 
x.mesh loc=20 spac=2.5 
x.mesh loc=250 spac=62.5 
#Y—-Mesh 
#Top contact 
































y.mesh loc=-211.1 spac=0.0125 

y-mesh loc=-211 spac=0.0125 

#Window 

y.mesh loc=-211 spac=0.0625 

y.mesh loc=-210.5 spac= 0.057293 

# Emitter 

y.mesh loc=-210.0417 spac=0.0052075 

# Emitter Depletion Width 

y.mesh loc=-210 spac=0.0052075 

# Base Depletion Width 

y-mesh loc=-209.5834 spac=0.052075 

# Base 

y.mesh loc=-205.2083 spac=1.1979 

y.mesh loc=-200 spac=1.1979 

# Substrate 

y.-mesh loc=0 spac=25 

#Regions 

# Contact 

region num=1 material=Gold x.min=-20 x.max=20 y.min=-211.1 \ 
y.max=-211 

region num=2 material=Vacuum x.min=-250 x.max=-20 y.min=-211.1 \ 
y.max=-211 

region num=3 material=Vacuum x.min=20 x.max=250 y.min=-211.1 \ 
y.max=-211 

# Window 

region num=4 material=AlGaAs x.comp=0.95 x.min=-250 x.max=250 \ 
y.min=-211 y.max=-210.5 

# Emmitter 
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region num=5 material=GaAs x.min=-250 x.max=250 y.min=-210.5 \ 
y.max=-210.0417 
# Emmitter Depletion Region 
region num=6 material=GaAs x.min=-250 x.max=250 y.min=-210.0417 \ 
y.max=-210 
# Base Depletion Region 
region num=7 material=GaAs x.min=-250 x.max=250 y.min=-210 \ 
y.max=-209.5834 
# Base 
region num=8 material=GaAs x.min=-250 x.max=250 y.min=—-209.5834 \ 
y.max=-200 
# Substrate 
region num=9 material=GaAs x.min=-250 x.max=250 y.min=-200 \ 
y.-max=0 
#Electrodes 

lectrode name=anode material=Gold x.min=-20 x.max=20 \ 
y.min=-211.1 y.max=-211 

lectrode name=cathode material=Gold x.min=-250 x.max=250 \ 
y.-min=0 y.max=0 

#Doping 

# Window 
doping uniform region=4 p.type conc=1000000000000000000 

# Emitter 
doping uniform region=5 p.type conc=1000000000000000000 
# Emitter Depletion Region 
doping uniform region=6 p.type conc=1000000000000000000 
# Base Depletion Region 
doping uniform region=7 n.type conc=100000000000000000 
# Base 
doping uniform region 
# Substrate 
doping uniform region=9 n.type conc=500000000000000000 
### Material properties ### 
## Vacuum 
material material=Vacuum real.index=3.3 imag.index=0 
## GaAs (from 
http://www.ioffe.rssi.ru/SVA/NSM/Semicond/GaAs/index.html1) 
material material=GaAs INDEX.FILE=opt/GaAsp.opt EG300=1.424 
#PERMITTIVITY=12.9 AFFINITY=4.07 NC300=4.7e17 NV300=9e18 \ 
COPT=7.2e-10 NSRHN=5e18 NSRHP=5e18 TAUNO=5e-9 TAUP0=3e-6 
## AlGaAs x=0.95 












































8 n.type conc=100000000000000000 

















































































































































































































































































































material material=AlGaAs INDEX.FILE=opt/AlGaAsx95.opt 
#EG300=1.5487 PERMI IVITY=12.616 AFFINITY=3.96 NC300=1.39e18 \ 
NV300=9.78e18 COPT=1.8e-10 TAUN0=4.5e-9 TAUP0=2e-8 

## Gold 
material material=Gold INDEX.FILE=opt/Gold.opt 

Ht tH ## Hitt Hitt Hitt # Het HHH HHH Hatt HEH HEH HE HH 

# Radiation effects # 

# Prior to Irradiation # 
Hitt Hee HEH HE Hit # HEE H Hate HEH HHH Hate HEH HE HEE HH 
trap acceptor e.level=0.71 density=3.0e11 degen=2 sign=4e-013 \ 
sigp=4e-015 region=5 

trap acceptor e.level=0.71 density=3.0el degen=2 sign=4e-013 \ 
sigp=4e-015 region=6 

# Models 

nodels region=4 ANALYTIC 

odels region=5 CONMOB BGN SRH 
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model 
model 
model 
model 


region=6 
region=7 
region=8 
region=9 


nnn nN DN 


models OPTR pri 
#Light beams 
beam num=1 x.origi 


rigin=-710.5 





power.file=AMOnrel. 
wavel.num=630 
beam num=2 x.origi 
power. file=AMOnrel 








-O 


wavel.start=0. 





rigin=-710.5 


wavel.start=0. 





wavel.num=825 
##Structure File 
struct outfile=st 
log off 

tonyplot 
HEHEHE HH HSH 
# Solves 
# 
## 














## 
and 


# 
Pp 


## 
lots 


# 











Hitt HH 
solve init 
set I=0.95 
method gummel newto 
solve b1=0.9 b2=0.9 


Ht tH 




















ructu 


structure/GaAs_ 
oun 
light 


Ht 


Ht 


Nn maxt 


re/ 


i 





# 





GaAs_L 





KL_O.str 
HEHEHE HHH HH 
IV curve 








Hit # tt 

















raps=10 it 


method newton maxtraps=10 itlimit=100 








ve b1=$I b2=$I 
ract name="isc" 
isc=Sisc*200000 
i1=Sisc/10 
12=$11+Sisc/10 
13=S$i2+Sisc/10 
14=S$i3+Sisc/10 
15=Si4+Sisc/10 
16=S$i5+Sisc/20 
17=Si6+Sisc/20 
18=Si7+Sisc/20 
19=Si8+Sisc/20 
i110=$i9+Sisc/20 
11=$i10+Sisc/40 
12=$1i11+$isc/40 
13=$112+$isc/40 
14=S$i13+Sisc/40 
15=$1i14+$isc/40 
16=$i15+Sisc/80 
17=S$i16+Sisc/80 
18=S$i1l7+Sisc/80 
119=$i118+Sisc/80 
120=$119+Sisc/80 
121=$i20+Sisc/80 
122=$i21+Sisc/80 
123=$122+Sisc/80 
set 124=$123+Sisc/80 
set i25=$124+Sisc/80 
##Luminous Curve 
contact name=anod 


sol 
ext 
set 
set 
set 
set 
set 
set 
set 
set 
set 
set 
set 
set 
set 
set 
set 
set 
set 
set 
set 
set 
set 
set 
set 
set 





ee oe 














max(i."cathode") 


current 


method newton maxtraps=10 itlimit=100 
log outfile=output/GaAs_LKL_0.log 
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Ht 


tt 


angle=90 back.refl \ 
12 wavel.end=0.75 \ 


angle=90 back.refl \ 
7501 wavel.end=2.4 \ 


Har aR aE aE aE aE aE aE aE aE aE aE EEE EEE 
# 
# 
Hear aE ae at aE aE aE aE aE aE aE aE aE EE HE 





limit=25 















































































































































solve ianode=-$i25 bl1=SI b2=SI 

solve ianode=-$i24 b1=SI b2=SI 

solve ianode=-$i23 b1=SI b2=SI 

solve ianode=-$i22 b1=SI b2=SI 

solve ianode=-$i21 b1=SI b2=SI 

solve ianode=-$i20 b1l=SI b2=SI 

solve ianode=-$il19 b1=SI b2=SI 

solve ianode=-$il18 bl1=SI b2=SI 

solve ianode=-$il17 b1=SI b2=SI 

solve ianode=-$il6 b1l=SI b2=SI 

solve ianode=-$i15 bl1=SI b2=SI 

solve ianode=-$il14 b1=SI b2=SI 

solve ianode=-$i13 b1=SI b2=SI 

solve ianode=-$il12 b1=SI b2=SI 

solve ianode=-$ill b1=SI b2=SI 

solve ianode=-$il0 b1l=SI b2=SI 

solve ianode=-$i9 b1=SI b2=SI 

solve ianode=-$i8 b1=SI b2=SI 

solve ianode=-$i7 b1=SI b2=SI 

solve ianode=-$i6 b1=SI b2=SI 

solve ianode=-$i5 bl1=SI b2=SI 

solve ianode=-$i4 b1=SI b2=SI 

solve ianode=-$i3 b1=SI b2=SI 

solve ianode=-$i2 b1=SI b2=SI 

solve ianode=-$il b1=SI b2=SI 

solve ianode=0 b1l=SI b2=SI 

log off 

tonyplot output/GaAs_LKL_0.log -set ivcurve.set 
HH Ht HH HH HH HH SE EH HH EE EH HH HE EE EE HH HE EE EH OH OE EEE EH OE HE EE 
# Solves and plots the dark IV curve 
# 

HHEEEEEEEE EE HH HH HE HF HH EH HE HE HE HE HE HE HE HH HE HE EE HE HH EH HH HE EH HE HH HH 
solve init 

method newton maxtraps=10 itlimit=100 CLIM.DD=1le8 
contact name=anode 

set Idrk= 

log outfile=output/drkGaAs_LKL_0.log 
solve vanode=0.00 bl=$Idrk b2=SIdrk 
solve vanode=0.20 bl=S$Idrk b2=SIdrk 
solve vanode=0.40 bl=$Idrk b2=SIdrk 
solve vanode=0.60 bl=$Idrk b2=SIdrk 
solve vanode=0.80 bl=SIdrk b2=STIdrk 
solve vanode=0.82 bl=$Idrk b2=SIdrk 
solve vanode=0.84 bl=SIdrk b2=SIdrk 
solve vanode=0.86 bl=SIdrk b2=STIdrk 
solve vanode=0.88 bl=$Idrk b2=SIdrk 
solve vanode=0.90 bl=SIdrk b2=STIdrk 
solve vanode=0.91 bl=S$Idrk b2=SIdrk 
solve vanode=0.92 bl=S$Idrk b2=SIdrk 
solve vanode=0.93 bl=SIdrk b2=SIdrk 
solve vanode=0.94 bl=S$Idrk b2=SIdrk 
solve vanode=0.95 bl=S$Idrk b2=SIdrk 
solve vanode=0.96 bl=S$Idrk b2=SIdrk 
solve vanode=0.97 bl=$Idrk b2=SIdrk 
solve vanode=0.98 bl=SIdrk b2=STIdrk 
solve vanode=0.99 bl=S$Idrk b2=SIdrk 
solve vanode=1.00 bl=S$Idrk b2=SIdrk 
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Hat HT Ht 
tt 
# 
Hat HT Ht 

















solve vanode=1.01 bl=$Idrk b2=SIdrk 

solve vanode=1.02 bl=S$Idrk b2=SIdrk 

solve vanode=1.04 bl=$Idrk b2=SIdrk 

solve vanode=1.06 bl=$Idrk b2=SIdrk 

solve vanode=1.08 bl=$Idrk b2=SIdrk 

solve vanode=1.10 bl=$Idrk b2=SIdrk 

solve vanode=1.13 bl=$Idrk b2=SIdrk 

solve vanode=1.16 bl=S$Idrk b2=SIdrk 

solve vanode=1.19 bl=S$Idrk b2=SIdrk 

solve vanode=1.22 b1l=S$Idrk b2=SIdrk 

solve vanode=1.26 bl=S$Idrk b2=SIdrk 

solve vanode=1.32 b1l=S$Idrk b2=SIdrk 

solve vanode=1.40 bl=$Idrk b2=SIdrk 

log off 

tonyplot output/drkGaAs_LKL_0.log 

HH HH Ht HE HH HH SE EH HH Ht EE HH HE EE HH HH 

# Frequency log provides output for 

# Structure files provide pictures 

# recombination of carriers within the structure. 
HH Ht HH HH HH HE EE HH HE EE EH HH EE EH HH EF 

set I1=1.0 

output con.band val.band 

solve init 

log outfile= output/GaAs_LKL_0_freq.log 
contact name=anode 

struc outfile=structure/illum/LKLO-0.20.str 
solve vanode=0 bl1=SI1 b2=SI1 lambda=0.2 
struc outfile=structure/illum/LKLO-0.22.str 
solve vanode=0 b1l=SI1 b2=S$I1 lambda=0.22 
struc outfile=structure/illum/LKLO-0.24.str 
solve vanode=0 b1l=SI1 b2=SI1 lambda=0.24 
struc outfile=structure/illum/LKLO-0.26.str 
solve vanode=0 b1l=SI1 b2=SI1 lambda=0.26 
struc outfile=structure/illum/LKLO-0.28.str 
solve vanode=0 b1=SI1 b2=SI1 lambda=0.28 
struc outfile=structure/illum/LKLO-0.30.str 
solve vanode=0 b1=SI1 b2=SI1 lambda=0.3 
struc outfile=structure/illum/LKLO-0.32.str 
solve vanode=0 b1l=SI1 b2=SI1 lambda=0.32 
struc outfile=structure/illum/LKLO-0.33.str 
solve vanode=0 b1=SI1 b2=SI1 lambda=0.33 
struc outfile=structure/illum/LKLO-0.335.str 
solve vanode=0 b1=SI1 b2=SI1 lambda=0.335 
struc outfile=structure/illum/LKLO-0.34.str 
solve vanode=0 b1l=SI1 b2=S$I1 lambda=0.34 
struc outfile=structure/illum/LKLO-0.35.str 
solve vanode=0 b1=SI1 b2=SI1 lambda=0.35 
struc outfile=structure/illum/LKLO-0.36.str 
solve vanode=0 b1=SI1 b2=SI1 lambda=0.36 
struc outfile=structure/illum/LKLO-0.38.str 
solve vanode=0 b1=SI1 b2=SI1 lambda=0.38 
struc outfile=structure/illum/LKLO-0.40.str 
solve vanode=0 bl1=SI1 b2=SI1 lambda=0.4 
struc outfile=structure/illum/LKLO-0.42.str 
solve vanode=0 b1l=SI1 b2=SI1 lambda=0.42 
struc outfile=structure/illum/LKLO-0.44.str 
solve vanode=0 b1l=SI1 b2=S$I1 lambda=0.44 
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curv 


Hat Ha aE HE HE HE HE EE EE EE 
spectral respons 
of photogeneration and 


tt 
# 
i 


HE aE aE aE HE aE EE EE EE EE EE EE EE 


structure/illum/LKLO-0.46.str 
1=$I1 b2=SI1 lambda=0.46 
ructure/illum/LKLO-0.48.str 
1=SI1 b2=SI1 lambda=0.48 
ructure/illum/LKLO-0.50.str 
1=SI1 b2=SI1 lambda=0.5 
ructure/illum/LKLO-0.52.str 
1=SI1 b2=SI1 lambda=0.52 
ructure/illum/LKLO-0.54.str 
1=SI1 b2=S$I1 lambda=0.54 
ructure/illum/LKLO-0.56.str 
1=SI1 b2=SI1 lambda=0.56 
ructure/illum/LKLO-0.58.str 
1=SI1 b2=S$I1 lambda=0.58 
ructure/illum/LKLO-0.60.str 
1=SI1 b2=SI1 lambda=0.6 
ructure/illum/LKLO-0.62.str 
1=SI1 b2=S$I1 lambda=0.62 
ructure/illum/LKLO-0.64.str 
1=SI1 b2=SI1 lambda=0.64 
ructure/illum/LKLO-0.66.str 
1=SI1 b2=SI1 lambda=0.66 
ructure/illum/LKLO-0.68.str 
1=SI1 b2=SI1 lambda=0.68 
ructure/illum/LKLO-0.70.str 
1=SI1 b2=S$I1 lambda=0.7 
ructure/illum/LKLO-0.72.str 
1=SI1 b2=S$I1 lambda=0.72 
ructure/illum/LKLO-0.74.str 
1=SI1 b2=S$I1 lambda=0.74 
ructure/illum/LKLO-0.76.str 
1=SI1 b2=S$I1 lambda=0.76 
struc ou ructure/illum/LKLO-0.78.str 
solve va 1=SI1 b2=SI1 lambda=0.78 
struc outfile=structure/illum/LKLO-0.80.str 
solve va 1=SI1 b2=SI1 lambda=0.8 
struc ou ructure/illum/LKLO-0.82.str 
solve va 1=SI1 b2=SI1 lambda=0.82 
struc ou ructure/illum/LKLO-0.84.str 
solve va 1=SI1 b2=S$I1 lambda=0.84 
struc ou ructure/illum/LKLO-0.86.str 
solve va 1=SI1 b2=SI1 lambda=0.86 
struc ou cucture/illum/LKLO-0.88.str 
solve va 1=SI1 b2=SI1 lambda=0.88 
struc outfile=structure/illum/LKLO-0.89.str 
solve vanode=0 b1l=SI1 b2=SI1 lambda=0.89 
struc out structure/illum/LKLO-0.90.str 
solve vanode b1=$I1 b2=SI1 lambda=0.9 

log off 


struc ou 
solve va 
struc ou 
solve va 
struc ou 
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solve va 
struc ou 
solve va 
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solve va 
struc ou 
solve va 
struc ou 
solve va 
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solve va 
struc ou 
solve va 
struc ou 
solve va 
struc ou 
solve va 
struc ou 
solve va 
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Pt sat sat set eat sat sate spate StS tet sad 
cose a |e 9 a | ee a | > | 


ode 
struc ou 
solve va 
struc ou 
solve va 
struc ou 
solve va 
struc ou 
solve va 





th 
H 
ll @ 
n 
ct 


ololotloto tl 
n 





i 
il o 





n 





e 
il o 


ode 


n 





es 
il 


ode 


n 
OF Ct sO sat Oeste ch AG 





P TYP TOT OT oO 


bs 
il © 





ll @O 


n 





co 
il o 


n 





2 
il o 


n 





e 
il o 


n 
Tt otc otc oct o 





patrsctsadcdsoadcsodc 
(e) 
Q 
oO 
ololhlotlototl 


ts 
il o 



























































2. TRAP STATEMENTS FOR SIMULATED FLUENCE LEVELS 
The following trap statements can be substituted within the input deck listed in 
Section 1. These traps represent the effects of 1 MeV electron radiation for fluence le vels 
of 10'4e/cm’, 10° e/cm’, 10!%e/cm’. 
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sign=2.6 


015 region=7 
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Hat Ht a EH HE aE aE HE aE EE EEE EEE EEE EEE EEE EEE EEE EE EEE EEE EE EEE EE EE EE 

# Radiation effects 

# 1 MeV electron Fluence level 10to014 e/cm2 

Hat oH HEH EH EE EE HEE EEE EEE EE EEE EE EE EE EE EE EE EE EE EE EE 
trap acceptor e.level=0.71 density=3.0e12 degen=2 sign=4e-013 \ 
sigp=4e-015 region=5 

trap acceptor e.level=0.71 density=3.0e12 degen=2 sign=4e-013 \ 
sigp=4e-015 region=6 

trap donor e.level=0.14 density=1.8e13 degen=2 sigp=5.7e-013 \ 
sign=5.7e-015 region=7 

trap donor e.level=0.41 density=8.2e13 degen=2 sigp=2.6e-013 \ 
sign=2.6e-015 region=7 

trap donor e.level=0.71 density=3.0e10 degen=2 sigp=8.3e-013 \ 
sign=8.3e-015 region=7 

trap donor e.level=0.90 density=8.8el11 degen=2 sigp=3e-011 \ 
sign=3e-013 region=7 

trap donor e.level=0.14 density=1.8e13 degen=2 sigp=5.7e-013 \ 
sign=5.7e-015 region=8 

trap donor e.level=0.41 density=8.2e13 degen=2 sigp=2.6e-013 \ 
sign=2.6e-015 region=8 

trap donor e.level=0.71 density=3.0e10 degen=2 sigp=8.3e-013 \ 
sign=8.3e-015 region=8 

trap donor e.level=0.90 density=8.8el11 degen=2 sigp=3e-011 \ 
sign=3e-013 region=8 

trap donor e.level=0.14 density=1.8e13 degen=2 sigp=5.7e-013 \ 
sign=5.7e-015 region=9 

trap donor e.level=0.41 density=8.2e13 degen=2 sigp=2.6e-013 \ 
sign=2.6e-015 region=9 

trap donor e.level=0.71 density=3.0e10 degen=2 sigp=8.3e-013 \ 
sign=8.3e-015 region=9 

trap donor e.level=0.90 density=8.8el11 degen=2 sigp=3e-011 \ 
sign=3e-013 region=9 

eat oH EH EH EH EH EE HE EE EE EE OE EE EEE EE EE EE EE EE EE EE 

# Radiation effects # 

# 1 MeV electron Fluence level 10tol5 e/cm2 # 
Hatt HEH HE HEE Hat HH EH HE HE HE HE HEE EE OE EE HE EE EE EE HEE 
trap acceptor e.level=0.13 density=2.2e14 degen=2 sign=2.6e-015 
Ssigp=2.6e-017 region=5 

trap acceptor e.level=0.29 density=4.0e14 degen=2 sign=2.6e-013 
sigp=2.6e-015 region=5 

trap acceptor vel=0.35 density=8.0e13 degen=2 sign=6.4e-015 
sigp=6.4e-017 region=5 

trap acceptor e.level=0.71 density=6.4e13 degen=2 sign=4e-013 \ 
sigp=4e-015 region=5 

trap acceptor level=0.13 density=2.2e14 degen=2 sign=2.6e-015 
sigp=2.6e-017 region=6 

trap acceptor vel=0.29 density=4.0e14 degen=2 sign=2.6e-013 
sigp=2.6e-015 region=6 

trap acceptor vel=0.35 density=8.0e13 degen=2 sign=6.4e-015 
sigp=6.4e-017 region=6 

trap acceptor e.level=0.71 density=6.4e13 degen=2 sign=4e-013 \ 
sigp=4e-015 region=6 

trap donor e.level=0.41 density=2.1e13 degen=2 sigp=2.6e-013 \ 


ra 
ig 
ra 
ig 
ra 
ig 
ra 
ig 
ra 
ig 
ra 
ig 
ra 
ig 
ra 
ig 
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p donor e.level=0.71 density=1.7e13 degen=2 sigp=8.3e-013 \ 
n=8.3e-015 region=7 
p donor e.level=0.90 density=2.8e13 degen=2 sigp=3e-011 \ 
n=3e-013 region=7 
p donor e.level=0.41 density=2.1le13 degen=2 sigp=2.6e-013 \ 
n=2.6e-015 region=8 
p donor e.level=0.71 density=1.7e13 degen=2 sigp=8.3e-013 \ 
n=8.3e-015 region=8 
p donor e.level=0.90 density=2.8e13 degen=2 sigp=3e-011 \ 
n=3e-013 region=8 
p donor e.level=0.41 density=2.1e13 degen=2 sigp=2.6e-013 \ 
n=2.6e-015 region=9 
p donor e.level=0.71 density=1.7e13 degen=2 sigp=8.3e-013 \ 
n=8.3e-015 region=9 
p donor e.level=0.90 density=2.8e13 degen=2 sigp=3e-011 \ 
n=3e-013 region=9 
HHH HHH Ht Ht ae aE HE EH EE EE EE HEE EE EE EE EE EE EE EE EE 
adiation effects # 
MeV electron Fluence level 10tol6 e/cm2 # 
HHH HE HE HE EH EE HEE HE EE EE EE EE EE EE EE EE EE EE HH 
p acceptor e.level=0.13 density=8.4e14 degen=2 sign=2.6e-015 
p=2.6e-017 region=5 
p acceptor level=0.29 density=1.6e15 degen=2 sign=2.6e-013 
p=2.6e-015 region=5 
p acceptor e.level=0.35 density=1.0e15 degen=2 sign=6.4e-015 
p=6.4e-017 region=5 
p acceptor e.level=0.71 density=2.7e14 degen=2 sign=4e-013 \ 
p=4e-015 region=5 
p acceptor level=0.13 density=8.4e14 degen=2 sign=2.6e-015 
p=2.6e-017 region=6 
p acceptor vel=0.29 density=1.6e15 degen=2 sign=2.6e-013 
p=2.6e-015 region=6 
p acceptor vel=0.35 density=1.0e15 degen=2 sign=6.4e-015 
p=6.4e-017 region=6 
p acceptor e.level=0.71 density=2.7e14 degen=2 sign=4e-013 \ 
p=4e-015 region=6 
p donor e.level=0.41 density=8.8e13 degen=2 sigp=2.6e-013 \ 
n=2.6e-015 region=7 
p donor e.level=0.71 density=5.0e13 degen=2 sigp=8.3e-013 \ 
n=8.3e-015 region=7 
p donor e.level=0.90 density=6.5e13 degen=2 sigp=3e-011 \ 
n=3e-013 region=7 
p donor e.level=0.41 density=8.8e13 degen=2 sigp=2.6e-013 \ 
n=2.6e-015 region=8 
p donor e.level=0.71 density=5.0e13 degen=2 sigp=8.3e-013 \ 
n=8.3e-015 region=8 
p donor e.level=0.90 density=6.5e13 degen=2 sigp=3e-011 \ 
n=3e-013 region=8 
p donor e.level=0.41 density=8.8e13 degen=2 sigp=2.6e-013 \ 
n=2.6e-015 region=9 
p donor e.level=0.71 density=5.0e13 degen=2 sigp=8.3e-013 \ 
n=8.3e-015 region=9 
p donor e.level=0.90 density=6.5e13 degen=2 sigp=3e-011 \ 
n=3e-013 region=9 
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APPENDIX B. MATLAB SCRIPTS AND FUNCTIONS 


MATLAB scripts and functions were used to produce input decks that could be 
modified and accurately reproduced during iterative solutions [Ref. 5]. During the re- 
search of this method of radiation modeling, it was import to have the capability to ma- 
nipulate parameters that would now appear to be standardized. As such, these codes are 
not expressed in the most efficient method one would use model solar cells on a more 
broad scale. Section 1 is the MATLAB script to produce a comprehensive input deck 
specifying desired fluence levels and output requirements. Section 2 is a function that 
produces the bulk of the standard ATLAS input deck. Section 3 contains the functions to 
produce the lines of code that represent the various fluence levels. 


1. MATLAB CODE FOR PRIMARY INPUT FILE 
The primary script file used to manipulate the desired output is LKLcelll.m. This 


file allows the designation of trap settings, fluence levels and output curves and data. 


format long 


fclose(‘all’); 

clear; 

cle; 

rdlvl='allsets'; % 0, 14, 15, 16 or ‘allsets' 

t="b'; % ‘a’ donors in the p region or 'b' acceptors in the p region 

capt='cec'; % default 'chc' capture holes in p region or 'cec' capture electrons in the p-region 
fac=le-2; % multiplication factor of capture cross section sigma 


IVcurve=1;  _% | produce IV output 
drkcurve=1; = % 1 produce dark current output 
specurve=1; % | produce spectral response output 
file=char(' '); 
TRAPMOD='SRH'; 
if rdlvl=='allsets' 
for r=1:4 
rd=[0 14 15 16]; 
radlvl=rd(r); 
if radlvl==0 
radpart=LKLradO(t,capt,fac); 
fileO=LKLSilv 1 (radlvl,radpart,t,capt,fac,.[Vcurve,drkcurve,specurve, TRAPMOD); 
file=char(file, fileO); 
r=r+1; 
end 
if radlvl==14 
radpart=LKLrad14(t,capt,fac); 
filel4=LKLSilv 1 (radlv1,radpart,t,capt,fac,[Vcurve,drkcurve,specurve, TRAPMOD); 
file=char(file, file14); 
r=r+1; 
end 
if radlvl==15 
radpart=LKLrad15(t,capt,fac); 
file 15=LKLSilv 1 (radlv1radpart,t,capt,fac,1Vcurve,drkcurve,specurve, TRAPMOD); 
file=char(file, file 15); 
r=r+1; 
end 
if radlvl==16 
radpart=LKLrad16(t,capt,fac); 
filel16=LKLSilv 1 (radlv1 radpart,t,capt,fac,.[Vcurve,drkcurve,specurve, TRAPMOD); 


a3 


file=char(file, file 16); 
end 
end 
end 


radlvl=rdlvl; 
if radlvl== 
radpart=LKLradO(t,capt,fac); 
file=LKLSilv 1 (radlvLradpart,t,capt,fac,.[Vcurve,drkcurve,specurve, TRAPMOD) 
end 
if radlvl==14 
radpart=LKLrad14(t,capt,fac); 
file=LKLSilv | (radlv1,radpart,t,capt,fac,1Vcurve,drkcurve,specurve, TRAPMOD) 
end 
if radlvl==15 
radpart=LKLrad 15(t,capt,fac); 
file=LKLSilv 1 (radlvl,radpart,t,capt,fac,1Vcurve,drkcurve,specurve, TRAPMOD) 
end 
if radlvl==16 
radpart=LKLrad 1 6(t,capt,fac); 
file=LKLSilv | (radlv1,radpart,t,capt,fac,.1Vcurve,drkcurve,specurve, TRAPMOD) 
end 





% overwrite previous input deck 
infile=char(['LKL_AII']); 
f=fopen([infile '.in'],'w'); 


[x,y]=size(file); 

for i=1:x 
fprintf(f,'%s\n',deblank(file(,:))); 

end 


fclose(f); 


2 MATLAB CODE TO PRODUCE THE BULK SILVACO CODE 
The bulk of the ATLAS code does not change with fluence levels. The solar cell 
structure, material parameters, and calculative methods are kept standardized across solar 


cell models. The LKLSi/v1.m function produces this portion of the input deck. 


function file=LKLSilvl (radlvl, radpart,t, capt, fac, IlVcurve, drkcurve, specurve, TRAPMOD) 
SPARAMETERS 

fclose('all'); 

% GaAs Constants (external to Silvaco) 

% emmissivity (F/m) 

epsR=12.9; 

epso=8.8542e-17; 

% intrinsic 
ni=1.79e6; 
STemperature 
T=300; 
e=1.60218e-19; 
k=1.3806e-23; 


S%Cell structure 
SCell Width (um) 
cellW=500; 


SCap/Contact (%) optimum 8% 
contWper=8; 


SContact 
contTH=0.1; 
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SWindow 
windowTH=0.5; 
winDop=1el18; 


SEmmiter 
EmTH=0.5; 
Edop=1e18; 


SBase 
BaseTH=10; 
Bdop=1e17; 


SSubstrate 
subTH=200; 
subDop=5el7; 


2 


% Depletion Layer 


Vo=k*T* log (Bdop*Edop/ ( (ni) *2))/e; 

WDep= ( (2*epsR*epso*Vo* (BdoptEdop) ) / (e*Bdop*Edop) ) *.5*10%6; 
BDepTH=WDep/ (1+Bdop/Edop) ; 

EDepTH=WDep/ (1+Edop/Bdop) ; 


SStandard number of Divisions 
divs=8; 


% Cell Structure calculations ### 


ole 


cellWDiv=cellW/divs; 
width3d=100e6/cellWw; 
cellWHalf=cellW/2; 
contW=0.01*contWper*cellWHalf; 
contWDiv=contW/divs; 


subLo=0; 
subHi=subLo-subTH; 
subDiv=subTH/divs; 


BLo=subHi; 
BHi=BLo-BaseTH+BDepTH; 
BDiv=(BaseTH-BDepTH) /divs; 


BMid=BLo- (BaseTH+BDepTH) /2; 


BDepLo=BHi; 
BDepHi=BDepLo -BDepTH; 
BDepDiv=BDepTH/divs; 


EmDepLo=BDepHi; 
EmDepHi=EmDepLo-EDepTH; 
EmDepDiv=EDepTH/divs; 


EmLo=EmDepHi; 
EmHi=EmLo-EmTH+EDepTH; 
EmDiv=(EmTH-EDepTH) /divs; 





windowLo=EmHi; 
windowHi=windowLo-windowTH; 
windowDiv=windowTH/divs; 


cont Lo=windowHi; 
contHi=contLo-contTH; 
contDiv=contTH/divs; 
light Y=EmHi-500; 
file=char('go atlas',... 


['mesh width=' num2str(width3d) ], 
['#X-Mesh'], 
[ 


SP) 


\Cell is built in two mirrored halves +/- fo zero axis ie. 


( 


++@@*% 


'x.mesh loc=-' num2str(cellWHalf) ' spac=' num2str(cellWDiv)], 


| 


AAQ@QHt+ 


) 







































































'"x.mesh loc=—-' num2str(contW) ' spac=' num2str(contWDiv) ], 

'x.mesh loc=' num2str(contW) ' spac=' num2str(contWDiv) ], 

'x.mesh loc=' num2str(cellWHalf) ' spac=' num2str(cellWDiv)], 

'#Y¥—Mesh'], ... 

'#Top contact'], 

'y.mesh loc=' num2str(contHi) ' spac=' num2str(contDiv)], 

"'y.mesh loc=' num2str(contLo) ' spac=' num2str(contDiv)], 

'#Window'], : 

'y.mesh loc=' num2str(windowHi) ' spac=' num2str(windowDiv) ], 

"'y.mesh loc=' num2str(windowLo) ' spac= ' num2str(EmDiv)], 

'# Emitter'], ... 

'y.mesh loc=' num2str(EmLo) ' spac=' num2str(EmDepDiv) J], 

'# Emitter Depletion Width'], 

'y.mesh loc=' num2str(EmDepLo) ' spac=' num2str(EmDepDiv) ], 

'# Base Depletion Width'], 

"y.mesh loc=' num2str(BDepLo) ' spac=' num2str(BDepDiv) ], 

'# Base'], ... 

"'y.mesh loc=' num2str(BMid) ' spac=' num2str(BDiv)], 

'y.mesh loc=' num2str(BLo) ' spac=' num2str(BDiv)], 

'# Substrate'], 

"'y.mesh loc=' num2str(subLo) ' spac=' num2str(subDiv)], 

'#Regions'], 

'e: Contact."-],) 

"region num=l1 material=Gold x.min=-' num2str(contW) ' x.max="' num2str(contW) ' 
y-min=' num2str(contHi) ' y.max=' num2str(contLo) ], . 

"region num=2 material=Vacuum x.min=-' num2str(cellWHalf) ' x.max=—' num2str(contW) 
' -y.min=' num2str(contHi) ' y.max=' num2str(contLo)], 

"region num=3 material=Vacuum x.min="' num2str(contW) ' x.max=' num2str(cellWHalf) ' 
y-min=' num2str(contHi) ' y.max=' num2str(contLo)], 

'# Window'], 

"region num=4 material=AlGaAs x.comp=0.95 x.min=—' num2str(cellWHalf) ' x.max=' 
num2str(cellWHalf) ' y.min=' num2str(windowHi) ' y.max=' num2str(windowLo) ], 

'# Emmitter'], 

"region num=5 material=GaAs x.min=—-' num2str(cellWHalf) ' x.max=' num2str(cellWHalf) 
'-y.min=' num2str(EmHi) ' y.max=' num2str(EmLo) J], 

'# Emmitter Depletion Region'], 

"region num=6 material=GaAs x.min=-' num2str(cellWHalf) ' x.max=' num2str(cellWHalf) 
' -y.min=' num2str(EmDepHi) ' y.max=' num2str(EmDepLo) J, 

'# Base Depletion Region'], 

"region num=7 material=GaAs x.min=-' num2str(cellWHalf) ' x.max=' num2str(cellWHalf) 
' y.min=' num2str(BDepHi) ' y.max=' num2str(BDepLo) ], 

'# Base'], 

"region num=8 material=GaAs x.min=—' num2str(cellWHalf) ' x.max=' num2str(cellWHalf) 
' y.min=' num2str(BHi) ' y.max=' num2str(BLo)], 

'# Substrate'], as 

"region num=9 material=GaAs x.min=-' num2str(cellWHalf) ' x.max="' num2str(cellWHalf) 
' y.min=' num2str(subHi) ' y.max=' num2str(subLo)], 

'#Electrodes'], 

"electrode name=anode material=Gold x.min=—' num2str(contW) ' x.max=' num2str(contW) 
" y.min=' num2str(contHi) ' y.max=' num2str(contLo)], . 

"electrode name=cathode material=Gold x.min=-' num2str(cellWHalf) ' x.max=' 
num2str(cellWHalf) ' y.min="' num2str(subLo) ' y.max=' num2str(subLo) ], 

'#Doping'], 

'# Window'], 

‘doping uniform region=4 p.type conc=' num2str(winDop) ], 

'# Emitter'], 

‘doping uniform region=5 p.type conc=' num2str(Edop)], 

'# Emitter Depletion Region'], ... 

‘doping uniform region=6 p.type conc=' num2str(Edop)], 

'# Base Depletion Region'], 

‘doping uniform region=7 n.type conc=' num2str(Bdop)], 

'"# Base'], 

‘doping uniform region=8 n.type conc=' num2str(Bdop)], 

'# Substrate'], 

"doping uniform region=9 n.type conc=' num2str(subDop) ], 

'## Material properties ##'], 

"# Vacuum'], 

"material material=Vacuum real.index=3.3 imag.index=0'], 

'# GaAs (from http://www.ioffe.rssi .ru/SVA/NSM/Semicond/GaAs/index.html)'], 

"material material=GaAs INDEX.FILE=opt/GaAsp.opt EG300=1.424 '], 3 

'#PERMITTIVITY=12.9 AFFINITY=4.07 NC300=4.7e17 NV300=9e18 COPT=7.2e-10 \'], 
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power. 


power. 











'#NSRHN=5e18 NSRHP=5e18 TAUNO=5e-9 TAUP0=3e-6'], 

'## AlGaAs x=0.95'J, 

‘material material=AlGaAs INDEX.FILE=opt/AlGaAsx95.opt '], 
'#EG300=1.5487 PERMITTIVITY=12.616 AFFINITY=3.96 NC300=1. 39018 x4 
'#NV300=9.78e18 COPT=1.8e-10 TAUNO=4.5e-9 TAUP0=2e-8 '], 

'# Gold'], 

"material material=Gold INDEX.FILE=opt/Gold.opt'], 

'##Radiation effects'], 

radpart J], 

'# Models'], 

"models region=4 ANALYTIC' ], 3s 

"models region=5 CONMOB BGN ' TRAPMOD], 

"models region=6 CONMOB BGN ' TRAPMOD], 

‘models region=7 CONMOB ' TRAPMOD], 
"models region=8 CONMOB ' TRAPMOD], 
‘models region=9 CONMOB ' TRAPMOD], 
"models OPTR print'], 
'#Light beams'], 


"beam num=1 x.origin=0 y.origin="' num2str(lightY) ' angle=90 back.refl 
file=AMOnrel.spec wavel.start=0.12 wavel.end=0.75 wavel.num=630'], 
"beam num=2 x.origin=0 y.origin=' num2str(lightY) ' angle=90 back.refl 





file=AMOnrel.spec wavel.start=0.7501 wavel.end=2.4 wavel.num=825'], 
"##Structure File'], 

"struct outfile=structure/GaAs_LKL_' num2str(radlvl) '.str'], 

‘log off'], 

'#tonyplot structure/GaAs_LKL_' num2str(radlvl) '.str']); 


if IVcurve== 


file=char(file,'solve init', 
"set I=0.95', 
‘method gummel newton maxtraps=10 itlimit=25', 
"solve b1=0.9 b2=0.9', 
"method newton maxtraps=10 itlimit=100', 
"solve b1=SI b2=SI', 
"extract name="isc" max(i."cathode")', 
['set isc=Sisc*' num2str(width3d) ], 
"set il=Sise/10',; 
"set 12=$ilt+$Sisc/10', 
"set 13=$i2+$Sisc/10', 
"set i4=$1i13+$Sisc/10', 
"set 15=$i4+$Sisc/10', 
"set i6=$i5+$Sisc/20', 
"set i7=$i6+Sisc/20', 
"set 18=Si7+$Sisc/20', 
"set i9=$i8+Sisc/20', 
"set i110=$i9+Sisc/20', 
"set i11=$i10+Sisc/40', 
"set 112=$111+$isc/40', 
"set i13=$112+Sisc/40', 
"set i14=$i13+Sisc/40', 
"set i15=$i114+Sisc/40', 
"set i16=$115+Sisc/80', 
"set i17=$i16+Sisc/80', 
"set i18=$117+Sisc/80', 
"set i19=$118+Sisc/80', 
"set i20=$119+Sisc/80', 
"set i121=$i20+$Sisc/80', 
"set i22=$i121+Sisc/80', 
"set i123=$122+Sisc/80', 
"set i24=$i23+Sisc/80', 
"set i125=$124+Sisc/80', 
'##Luminous Curve', 
"contact name=anode current', 
"method newton maxtraps=10 itlimit=100', 
['log outfile=output/GaAs_LKL_' num2str(radlvl) '.log'], 
"solve ianode=-$i25 b1=SI b2=SI', 
"solve ianode=-$i24 b1l=SI b2=SI', 
"solve ianode=-$i23 b1=SI b2=SI', 
"solve ianode=-$i22 b1l=SI b2=SI', 
"solve ianode=-$i21 bl=SI b2=SI', 
"solve ianode=-$i20 b1=SI b2=SI', 


af 








"solve ianode=-$il19 
"solve ianode=-$i18 
"solve ianode=-$il7 
"solve ianode=-$il6 =SI b2=SI', 
"solve ianode=-$i15 =SI b2=SI', 


b1=SI b2=SI', 
b 
b 
b 
b 

"solve ianode=-$i14 b1=SI b2=SI', 
b 
b 
b 
b 


=SI b2=SI', 
=SI b2=SI', 


"solve ianode=-$i13 =f] b2=ST", 
"solve ianode=-$il12 =S$I b2=SI', 
"solve ianode=-$ill b1=$I b2=SI', 
"solve ianode=-$i10 =SI b2=SI', 
"solve ianode=-$i9 =SI b2=SI', 
"solve ianode=-$i8 =SI b2=SI', 
"solve ianode=-$i7 =SI b2=SI', 
"solve ianode=-$i6 =SI b2=SI', 














TOO COOOOOS 
Il 
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"solve ianode=-$i5 b2=SI', 
"solve ianode=-$i4 =SI b2=SI', 
"solve ianode=-$i3 =SI b2=SI', 
"solve ianode=-$i2 =SI b2=SI', 
"solve ianode=$il =SI b2=SI', 








"solve ianode=0 b1=$I b2=SI', 

‘log off', ‘ 

['#tonyplot output/GaAs_LKL_' num2str(radlvl) '.log -set ivcurve.set']); 
end 


if drkcurve== 

drkplot=1 

file=char(file,'solve init', ... 
"method newton maxtraps=10 itlimit=100 CLIM.DD=1e8', 
'##Solving Dark Current', 
"contact name=anode', 
"set Idrk=0', 
['log outfile=output/drkGaAs_LKL_' num2str(radlvl) '.log'], 

















"solve vanode=0.00 bl=SIdrk b2=SIdrk', 
"solve vanode=0.20 bl=SIdrk b2=SIdrk', 
"solve vanode=0.40 bl=SIdrk b2=SIdrk', 
"solve vanode=0.60 bl=SIdrk b2=SIdrk', 
"solve vanode=0.80 bl=SIdrk b2=SIdrk', 
"solve vanode=0.82 b1l=SIdrk b2=SIdrk', 
"solve vanode=0.84 bl=SIdrk b2=SIdrk', 
"solve vanode=0.86 bl=$Idrk b2=SIdrk', 
"solve vanode=0.88 bl=SIdrk b2=SIdrk', 
"solve vanode=0.90 bl=SIdrk b2=SIdrk', 
"solve vanode=0.91 bl=SIdrk b2=SIdrk', 
"solve vanode=0.92 b1l=SIdrk b2=SIdrk', 
"solve vanode=0.93 bl=S$Idrk b2=SIdrk', 
"solve vanode=0.94 bl=SIdrk b2=SIdrk', 
"solve vanode=0.95 bl=SIdrk b2=SIdrk', 
"solve vanode=0.96 bl=SIdrk b2=SIdrk', 
"solve vanode=0.97 bl=SIdrk b2=SIdrk', 
"solve vanode=0.98 bl=SIdrk b2=SIdrk', 
"solve vanode=0.99 bl=SIdrk b2=SIdrk', 
"solve vanode=1.00 bl=SIdrk b2=SIdrk', 
"solve vanode=1.01 bl=SIdrk b2=SIdrk', 
"solve vanode=1.02 bl=SIdrk b2=SIdrk', 
"solve vanode=1.04 bl=SIdrk b2=SIdrk', 
"solve vanode=1.06 bl=SIdrk b2=SIdrk', 
"solve vanode=1.08 bl=SIdrk b2=SIdrk', 
"solve vanode=1.10 bl=SIdrk b2=SIdrk', 
"solve vanode=1.13 bl=SIdrk b2=SIdrk', 
"solve vanode=1.16 bl=$Idrk b2=SIdrk', 
"solve vanode=1.19 bl=SIdrk b2=SIdrk', 
"solve vanode=1.22 bl=SIdrk b2=SIdrk', 
"solve vanode=1.26 bl=SIdrk b2=SIdrk', 
"solve vanode=1.32 b1l=SIdrk b2=SIdrk', 
"solve vanode=1.40 bl=S$Idrk b2=SIdrk', 
V Log of" y 

['tonyplot output/drkGaAs_LKL_' num2str(radlvl) '.log']); 





end 


if specurve== 
specplot=1 
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file=char (file, 'set I1=1.0', 
‘output con.band val.band', 





solve init 


"log outfil 
"contact name=anode', 
file=structure/ill 
b2=SI1 


"struc out 


solve vanode=0 b1=ST1 
file=structure/ill 


TStrue out 


solve vanode=0 b1=$ b2=S 
file=structure/ill 


"struc out 


solve vanode=0 b1=SI1 b2=SI1 
file=structure/ill 


‘struc: out 


solve vanode=0 b1=$ b2=S 
file=structure/ill 


*Struc out 


solve vanode=0 bl1=$SI1 b2=SI 
file=structure/ill 


"struc out 


solve vanode=0 b1=S] b2=SI 
file=structure/ill 


PSPS Ore 


solve vanode=0 b1=$I1 b2=$ 
file=structure/ill 
solve vanode=0 b1=STI1 b2=SI 
file=structure/ill 


"struc out 


"struc out 


solve vanode=0 b1=SI1 b2=ST1 
file=structure/ill 


"struc out 


solve vanode=0 b1=$ b2=SI 
file=structure/ill 
solve vanode=0 b1=$ b2=$ 
file=structure/ill 
solve vanode=0 b1=$I1 b2=$ 
file=structure/ill 


"struc out 


*SELUC: Out 


‘spruce out 


solve vanode=0 b1=SI1 b2=S$I 
file=structure/ill 


"struc out 


solve vanode=0 b1=S] b2=ST1 
file=structure/illum 
solve vanode=0 b1=$I1 b2=$ 
file=structure/ill 


"struc out 


"struc out 


solve vanode=0 b1i=$ b2=S J 2 
file=structure/illum/LKL' num2str(radlvl) 
solve vanode=0 b1=SI1 b2=SI a 
file=structure/ill 
solve vanode=0 b1=$ b2=$ 
file=structure/ill 
solve vanode=0 b1=SI1 b2=SI 
file=structure/ill 
solve vanode=0 b1=SI1 b2=SI1 
file=structure/ill 
solve vanode=0 b1=$ b2=$ 
file=structure/ill 
solve vanode=0 b1=S$I1 b2=S1 
file=structure/ill 
solve vanode=0 b1=$ b2=$ 
file=structure/ill 
solve vanode=0 b1=SI1 b2=SI 
file=structure/ill 
solve vanode=0 b1=$I1 b2=$ 
file=structure/ill 
solve vanode=0 b1=$ b2=SI 
file=structure/ill 
solve vanode=0 b1=SI1 b2=SI 
file=structure/ill 
solve vanode=0 b1=$ b2=$ 
file=structure/ill 
solve vanode=0 b1=$I1 b2=$T1 


Struc owe 


"struc out 


"struc out 


"etre out 


"struc out 


"struc out 


"struc out 


Tstru® our 


rstruc. out 


TStruc: out 


"struc out 


"struc out 


"struc out 


"struc out 


“struc out 


"struc out 





"struc out 





' 
rose. 






































file=structure/ill 
solve vanode=0 b1=$I1 b2=$ 
file=structure/ill 
solve vanode=0 b1=$ b2=$ 


le= output/GaAs_LKL_' num2str(radlvl) 








file=structure/illum 
solve vanode=0 b1=$ b2=$ 
file=structure/ill 
solve vanode=0 b1=$ b2=$ 
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'_freq.log'], 
um/LKL' num2str(radlvl) '-0.20.str' 
lambda=0.2', ... 
um/LKL' num2str(radlvl) '-0.22.str' 
lambda=0.22', ... 
um/LKL' num2str(radlvl) '-0.24.str' 
lambda=0.24', ... 
um/LKL' num2str(radlvl) '-0.26.str' 
lambda=0.26', ... 
um/LKL' num2str(radlvl) '-0.28.str' 
1 lambda=0.28', ... 
um/LKL' num2str(radlvl) '-0.30.str' 
lambda=0.3', ... 
um/LKL' num2str(radlvl) '-0.32.str' 
lambda=0.32', ... 
um/LKL' num2str(radlvl) '-0.33.str' 
lambda=0.33', ... 
um/LKL' num2str(radlvl) '-0.335.str 
lambda=0.335', ... 
um/LKL' num2str(radlvl) '-0.34.str' 
lambda=0.34', ... 
lum/LKL' num2str(radlvl) '-0.35.str' 
lambda=0.35', ... 
um/LKL' num2str(radlvl) '-0.36.str' 
lambda=0.36', ... 
um/LKL' num2str(radlvl) '-0.38.str' 
lambda=0.38', ... 
um/LKL' num2str(radlvl) '-0.40.str' 
lambda=0.4', ... 
/LKL' num2str(radlvl) '-0.42.str' 
lambda=0.42"', ... 
um/LKL' num2str(radlvl) '-0.44.str' 
lambda=0.44', . 
'"-0.46.str' 
lambda=0.46', ... 
um/LKL' num2str(radlvl) '-0.48.str' 
lambda=0.48', ... 
um/LKL' num2str(radlvl) '-0.50.str' 
lambda=0.5', ... 
um/LKL' num2str(radlvl) '-0.52.str' 
lambda=0.52', ... 
lum/LKL' num2str(radlvl) '-0.54.str' 
lambda=0.54', ... 
um/LKL' num2str(radlvl) '-0.56.str' 
lambda=0.56', ... 
um/LKL' num2str(radlvl) '-0.58.str' 
lambda=0.58', ... 
um/LKL' num2str(radlvl) '-0.60.str' 
1 lambda=0.6', ... 
um/LKL' num2str(radlvl) '-0.62.str' 
lambda=0.62', ... 
um/LKL' num2str(radlvl) '-0.64.str' 
1 lambda=0.64', ... 
um/LKL' num2str(radlvl) '-0.66.str' 
lambda=0.66', ... 
um/LKL' num2str(radlvl) '-0.68.str' 
lambda=0.68', ... 
um/LKL' num2str(radlvl) '-0.70.str' 
lambda=0.7', ... 
/LKL' num2str(radlvl) '-0.72.str' 
lambda=0.72', ... 
um/LKL' num2str(radlvl) '-0.74.str' 
lambda=0.74', ... 
um/LKL' num2str(radlvl) '-0.76.str' 
lambda=0.76', ... 
um/LKL' num2str(radlvl) '-0.78.str' 
lambda=0.78', 





























"struc outfile=structure/illum/LKL' num2str(radlvl) '-0.80.str'], 
"solve vanode=0 b1=$] b2=SI1 lambda=0.8', ... 

"struc outfile=structure/illum/LKL' num2str(radlvl) '-0.82.str'], 
"solve vanode=0 b1=S b2=$ lambda=0.82', ... 

"struc outfile=structure/illum/LKL' num2str(radlvl) '-0.84.str'], 
"solve vanode=0 bl=$ b2=$ lambda=0.84', 

"struc outfile=structure/illum/LKL' num2str(radlvl) '-0.86.str'], 
"solve vanode=0 b1=S b2=$ lambda=0.86', ... 

"struc outfile=structure/illum/LKL' num2str(radlvl) '-0.88.str'], 
"solve vanode=0 b1=S b2=$ lambda=0.88', ... 

"struc outfile=structure/illum/LKL' num2str(radlvl) '-0.89.str'], 
"solve vanode=0 b1=$ b2=$I1 lambda=0.89', ... 

"struc outfile=structure/illum/LKL' num2str(radlvl) '-0.90.str'], 
"solve vanode=0 b1=$] b2=SI1 lambda=0.9', 























‘log off'); 
end 


A? MATLAB CODE FOR RADIATION LEVELS 

For each level of radiation fluence a different set of trap statements must be used 
to represent the level of damage incurred. These functions are called by LKLcel/l.m ac- 
cording to the specified fluence level. 

a. Function LKLrad0.m 

The function LKLrad0.m produces the portion of the code that represents the 


grown-in defects present in the pre-irradiated solar cell model. 
function radpart=LKLrad0 (t,capt, fac) 


if t=='a'; 
typel='acceptor'; 
type2='donor'; 
elseif t=='b'; 
typel='donor'; 
type2='acceptor'; 
end 


xan? 
y='p'; 


if capt=='cec' 
x='"p'; 
y='n'; 

end 


sigl3=5.7e-13; 
sig41=2.6e-13; 
siga71=8.3e-13; 
sig90=3.0e-11; 
sigb71=4.0e-13; 
sigal3=num2str(fac*sigl3); 
siga4l=num2str(fac*sig41); 
sigaa7l=num2str(fac*siga71) ; 
siga90=num2str (fac*sig90); 
sigba71l=num2str(fac*sigb71) ; 


n=5; 


radpart=char(['trap ' type2 ' e.level=0.71 density=3.0e11 degen=2 sig' y '=' 
num2Zstr(sigb71) ' sig' x '=' sigba7l ' region=' num2str(n)]); 


n=6; 
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radpart=char(radpart,['trap ' type2 ' e.level=0.71 density=3ell degen=2 sig' y '=' 
num2Zstr(sigb71) ' sig' x '=' sigba71l ' region=' num2str(n)]); 


b. Function LKLrad14.m 
The function LKLrad14.m produces the portion of the code that represents the de- 


fects introduced by exposure to 1MeV electron radiation at a fluence level of 10!4e/cm”. 
function radpart=LKLrad14 (t,capt, fac) 


if t=='a'; 
typel='acceptor'; 
type2='donor'; 

elseif t=='b'; 
typel='donor'; 
type2='acceptor'; 


end 
sans 
y='p'; 


if capt=='cec' 
x='pl; 
y='n'; 

end 


sigl3=5.7e-13; 
sig41=2.6e-13; 
siga71=8.3e-13; 
sig90=3.0e-11; 
sigb71=4.0e-13; 


sigal3=num2str (fac*sigl3); 
siga4l=num2str(fac*sig41); 
sigaa7l=num2str(fac*siga71) ; 
siga90=num2str (fac*sig90); 
sigba71l=num2str (fac*sigb71) ; 


n=5; 

radpart=char(['trap ‘ type2 : e.level=0.71 density=3.0e12 degen=2 Sig" y '=! 
num2str(sigb71) ' sig! x '=' sigba71 ' region=' num2str(n)]); 

n=6; 

radpart=char(radpart,['trap ' type2 ' e.level=0.71 density=3.0e12 degen=2 sig' y '='! 
num2Zstr(sigb71) ' sig' x '=' sigba71l ' region=' num2str(n)]); 

for n=7:9 

radpart=char(radpart,['trap ' typel ' e.level=0.14 density=1.8e13 degen=2 sig' x '='! 
num2Zstr(sigl3) ' sig' y '=' sigal3 ' region=' num2str(n)], 

['trap ' typel ' e.level=0.41 density=8.2e13 degen=2 sig' x '=' num2str(sig41) ' sig' y 
"=" siga4l ' region=' num2str(n)]J, 

['trap ' typel ' e.level=0.71 density=3.0e10 degen=2 sig' x '=' num2str(siga71) ' sig' y 
"=" sigaa71 ' region=' num2str(n)], 

['trap ' typel ' e.level=0.90 density=8.8e1l1 degen=2 sig' x '=' num2str(sig90) ' sig' y 
"=" siga90 ' region="' num2str(n)]); 

end 


Cc. Function LKLrad15.m 
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The function LKLrad15.m produces the portion of the code that represents the de- 


fects introduced by exposure to 1MeV electron radiation at a fluence level of 10'*e/cm”. 
function radpart=LKLrad15(t, capt, fac) 

if t=="'a 
typel='acceptor'; 
type2='donor'; 
elseif t=='b'; 


typel='donor'; 

type2='acceptor'; 
end 
=a 
y='p'; 
if capt=='cec!' 

x='pl; 

at 1 

y= nj; 

end 


sig41=2.6e-13; 
siga71=8.3e-13; 
sig90=3.0e-11; 
sigl3=2.6e-15; 
sig29=2.6e-13; 
sig35=6.4e-15; 
sigb71=4.0e-13; 


siga4l=num2str(fac*sig41); 
sigaa7l=num2str(fac*siga71); 


siga90=num2s 


sigal3=num2st 


siga29=num2s 


siga35=num2st 


tr (fac*sig90); 
r(fac*sigl3); 
tr (fac*sig29); 
r(fac*sig35) 





sigba71l=num2str (fac*sighb71) ; 


n=5; 

radpart=char(['trap ' type2 ' e.level=0.13 density=2.2e14 degen=2 sig' y '=' 
num2Zstr(sigl3) ' sig' x '=' sigal3 ' region="' num2str(n)], 

['trap ' type2 ' e.level=0.29 density=4.0e14 degen=2 sig' y '=' num2str(sig29) ' sig' x 
"=" siga29 ' region=' num2str(n)]J, 

['trap ' type2 ' e.level=0.35 density=8.0e13 degen=2 sig' y '=' num2str(sig35) ' sig' x 
"=" siga35 ' region=' num2str(n)J, . 

['trap ' type2 ' e.level=0.71 density= 6.4e13 degen=2 sig' y '=' num2str(sigb71) ' sig' x 
'='  sigba71l ' region=' num2str(n)]); 

n=6; 

vradpart=char(radpart, ['trap ' type2 ' e.level=0.13 density=2.2e14 degen=2 sig' y '=' 
num2Zstr(sigl3) ' sig' x '=' sigal3 ' region=' num2str(n)], 

['trap ' type2 ' e.level=0.29 density=4.0e14 degen=2 sig' y '="' num2str(sig29) ' sig' x 
'='" siga29 ' region=' num2str(n)], 

['trap ' type2 ' e.level=0.35 denette 8.0e13 degen=2 sig' y '=' num2str(sig35) ' sig! x 
"=" siga35 ' region=' num2str(n)]J, 

['trap ' type2 ' e.level=0.71 densieye 6.4e13 degen=2 sig' y '=' num2str(sigb71) ' sig' x 
'='  sigba71l ' region=' num2str(n)]); 

for n=7:9 

radpart=char(radpart,['trap ' typel ' e.level=0.41 density=2.1le13 degen=2 sig' x '=' 
num2Zstr(sig41) ' sig' y '=' siga4l ' region=' num2str(n)], 

['trap ' typel ' e. Teves 0.71 density=1.7e13 degen=2 sig' x '=' num2str(siga71) ' sig' y 
"=" sigaa71 ' region=' num2str(n)], 

['trap ' typel ' e.level=0.90 density= 2. 8e13 degen=2 sig' x '=' num2str(sig90) ' sig' y 
'='  siga90 ' region=' num2str(n)]); 
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end 


d. Function LKLrad16.m 
The function LKLrad16.m produces the portion of the code that represents the de- 


fects introduced by exposure to 1MeV electron radiation at a fluence level of 10'*e/cm”. 
function radpart=LKLrad16(t,capt, fac) 


if t=="'a' 
typel='acceptor'; 
type2='donor'; 
elseif t==' 
typel='"donor'; 
type2='acceptor'; 
end 


ans 
y='p'i 


if capt=='cec' 
x='pl; 
yo'n'; 

end 


sig41=2.6e-13; 

siga71=8.3e-13; 
sig90=3.0e-11; 
sigl3=2.6e-15; 
sig29=2.6e-13; 
sig35=6.4e-15; 
sigb71=4.0e-13; 


siga4l=num2str (fac*sig41); 
sigaa7l=num2str(fac*siga71) ; 
siga90=num2str(fac*sig90) ; 
sigal3=num2str(fac*sig13) ; 
siga29=num2str (fac*sig29) ; 
siga35=num2str(fac*sig35) ; 
sigba7l=num2str (fac*sigb71) ; 





n=5; 


radpart=char(['trap ' type2 ' e.level=0.13 density=8.4e14 degen=2 sig' y '=' 


num2str(sig1l3) ' sig' x '=' sigal3 ' region="' num2str(n)], ' 

['trap ' type2 ' e.level=0.29 density=1.6e15 degen=2 sig' y '=' num2str(sig29) ' sig' x 
"=" siga29 ' region=' num2str(n)]J, 

['trap ' type2 ' e.level=0.35 densiiry= 1.0e15 degen=2 sig' y '=' num2str(sig35) ' sig' x 
'='" siga35 ' region=' num2str(n)], 

['trap ' type2 ' e.level=0.71 density=2.7e14 degen=2 sig' y '=' num2str(sigb71) ' sig' x 


'='  sigba71 ' region=' num2str(n)]); 
n=6; 


radpart=char(radpart, ['trap ' type2 ' e.level=0.13 density=8.4e14 degen=2 sig' y '=' 


num2str(sigl3) ' sig' x '=' sigal3 ' region=' num2str(n)], ‘ 

['trap ' type2 ' e.level=0.29 density=1.6e15 degen=2 sig' y '=' num2str(sig29) ' sig' x 

"=" siga29 ' region=' num2str(n)]J, F 

['trap ' type2 ' e.level=0.35 acneiy= 1.0e15 degen=2 sig' y '=' num2str(sig35) ' sig' x 

'=! eee " region=' num2str(n)], 4 

[ee " type2 ' e.level=0.71 densien= 2.7e14 degen=2 sig' y '=' num2str(sigb71) ' sig' x 
ahi 


Ser ' region=' num2str(n)]); 


for n=7:9 


radpart=char(radpart, ['trap ' typel ' e.level=0.41 density=8.8e13 degen=2 sig' x '=' 
num2str(sig41) ' sig' y '=' siga4l ' region=' num2str(n)], 
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[! 


rap ' typel ' e.level=0.71 density=5.0e13 degen=2 sig' 


sigaa7l 


rap ' typel ' e.level=0.90 density=6.5e13 degen=2 sig' x 


siga90 


region=' num2str(n)], 


region=' num2str(n)]); 
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x 


num2Zstr(siga71) 


num2str(sig90) 


sig' y 


sig' y 
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